Mon. Not. R. Astron. Soc. OOP. [Tlf38l C2011') Printed 6 October 2011 (MN WT^ style file v2.2) 



Testing pre-main sequence models: the power of a Bayesian 
approach. 

M. Gennaro^*t, P.G. Prada Moroni^'^, E.Tognelli^'^ 

^ Max-Planck-Institut fiir Astronomie, Kdnigstuhl 17, D-69117, Heidelberg, Germany 

^ Universitd di Pisa, Dipartimento di Fisica "E. Fermi", Largo B. Pontecorvo, 3, 1-56127, Pisa, Italy 

^Istituto Nazionale di Fisica Nucleare - Sezione di Pisa, Largo B. Pontecorvo, 3, 1-56127, Pisa, Italy 

Submitted to MNRAS 1848 February 21. Accepted 1945 April 25. 

ABSTRACT 

Pre-main sequence (PMS) models provide invaluable tools for the study of star forming 
regions as they allow to assign masses and ages to young stars. Thus it is of primary 
importance to test the models against observations of PMS stars with dynamically 
determined mass. We developed a Bayesian method for testing the present generation 
of PMS models which allows for a quantitative comparison with observations, largely 
superseding the widely used isochrones and tracks qualitative superposition. 

Using the available PMS data we tested the newest PISA PMS models establishing 
their good agreement with the observations. The data cover a mass range from 0.3 
to ^ S.IMq, temperatures from ^ 3 x 10'^ to ^ 1.2 x 10''' K and luminosities from 
~ 3 X 10~^ to ~ 60Lq. Masses are correctly predicted within 20% of the observed 
values in most of the cases and for some of them the difference is as small as 5%. 
Nevertheless some discrepancies are also observed and critically discussed. 

By means of simulations, using typical observational errors, we evaluated the 
spread of log Tsim — log Tree, i-©- simulated — recovered ages distribution of the single 
objects. We also found that stars in binary systems simulated as coeval might be 
recovered as non coeval, due to observational errors. The actual fraction of fake non 
coevality is a complex function of the simulated ages, masses and mass ratios. We 
demonstrated that it is possible to recover the systems' ages with better precision 
than for single stars using the composite age-probability distribution, i.e. the product 
of the components' age distributions. Using this valuable tool we estimated the ages 
of the presently observed PMS binary systems. 

Key words: Methods: statistical - binaries: general - Stars: fundamental parameters 
- Stars: pre-main sequence - binaries: eclipsing 



1 INTRODUCTION 

The current understanding of star formation processes 
largely relies on the ability of assigning ages and masses 
to young stars using pre-main sequence (PMS) models. The 
observed luminosity and effective temperature of stars in 
their early evolutionary stages can be translated into mass 
and age only by the comparison with PMS stellar tracks. 
Unfortunately the early evolutionary stages of the stellar 
life are among those less tightly constrained by observations 
and most uncertain from the theoretical point of view. This 
situation becomes progressively worse for stellar mass be- 
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low « 1.2 Mq. This is mainly a consequence of the poor 
treatment of superadiabatic convection. Moreover there are 
still large uncertainties on the main input physics describing 
the cold and dense matter typical of low-mass stars interiors 
adopted in modern evolutionary codes. This theoretical un- 
certainty is testified by the large discrepancy still present 
between different sets of low-mass PMS models (see e .g. 
ISiess et al.ll200d : iBaraffe et al.ll2002l : lTognem et al]|201ll ). 

An ever growing amount of detailed information is 
becoming availa ble for both sta r-forming regions in the 
Milky Way fsee [Reipurthll2q08allb | ) and in the Mag ellanic 
Clouds (|Cignoniet al.ll2009l . l20ialGouliermis et al.lboiol '). 
prompted by the remarkable improvement in the observa- 
tional techniques over the last decade. The aforementioned 
theoretical uncertainties imply that many of the properties 
inferred for these regions, such as the initial mass func- 
tion and the star formation history, strongly depend on 
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the adopted PMS models, particularly for stars less massive 
than 1.2 Mq. 

The importance of these kind of studies urges an empir- 
ical calibration of PMS tracks and isochrones based on a sta- 
tistically significant sample of young stars with precisely de- 
termined parameters (mass, temperature, radius, luminosity 
and chemical abundances). Most useful in this respect are 
the detached, double lined, eclipsing systems which directly 
provide stellar masses, temperatures and radii, yielding also 
a distance independent luminosity through the Stephan law. 
Unfortunately, the currently available sample of PMS stars 
in eclipsing binaries amounts only to 10 objects in 6 sys- 
tems (see Sect. [5] and Table |4]). Other important observa- 
tional constraints to PMS models are provided by astro- 
metric measurement of binary systems that can be resolved 
thanks to interferometric observations. For these systems the 
radii can not be determined though. Three astrometric bi- 
nary systems with both stars in the PMS phase are currently 
known and studied (Table [4]). A third technique providing 
mass values for young stars is based on spectroscopic obser- 
vations of circumstellar disks. From the keplerian velocities 
masses are inferred, although these measurements require 
an independent estimate of the distance to determine the 
linear value of the radius at which velocities are measured. 
Also the sample of known objects in this category is quite 
small with only 9 stars (see again Table |4]) 

As already shown in the early attempts to test PMS 
stellar t racks against observ a tions of eclipsing binaries 



(see e.g. IPalla fc StahW200ll; 'Hillenbrand fc White 



Stassun et al. 2004; Alccian et al. 2007; Bodcnetal 



2004 



2007 



Mathieu et al.ll2007 ). the mass values inferred from theoret- 



ical models are in reasonable agreement with the dynamical 
ones for intermediate mass stars, whereas for low-mass stars 
theoretical values tend to underestimate the stellar mass (see 
Fig. 3 in Mathieu ct al. 2007). 

From these studies it is also clear that the usefulness 
of such tests in constraining the theoretical PMS models is 
severely limited by the still scarce accuracy of the current 
empirical measurements of the other stellar parameters, i.e. 
the luminosity, the chemical a bundances and above all the 
effective temperatur e (see e.g. iHillenbrand fc White! [20041 : 
iMathieu et~al] 120071 ). In the next future both the size and 
the quality of the observed sample of test objects to calibrate 
PMS stellar tracks and isochrones are bound to increase. 

In the present paper we apply an objective Bayesian 
method to compare theoretical predictions with observa- 
tions, obtaining robust uncertainties for the output values 
and assessing the overall quality of the comparison. Since 
the method, which is detailed in Sect.[2l allows for the use 
of stellar tracks for a large and very fine grid of metallic- 
ity, mass and age values, we tested only our own Pisa-PMS 
models. These are calculated usin g the newest version of the 
FRANEC evolutionary code (see e.g. lTognelh et al.|[201ll '). The 
main characteristics of the models are described in Sect.O 
In Sect. |4] we asses the ability of the method to retrieve 
the stellar properties by means of synthetic tests. In Sect. [5] 
we describe the observational data set, which includes all 
the currently available low-mass PMS data. The complete 
data set is analysed in Sect.|6] where theoretical masses de- 
rived from our standard set of models are compared to the 
dynamical measurements. In Sect. [7] we compare the results 
for multiple sets of models. Section[8]is dedicated to detailed 



study of each binary system while the stars in the Taurus- 
Auriga association are analysed in Sect.|9] A summary with 
concluding remarks is presented in Sect. 1101 



2 THE BAYESIAN METHOD 

The general question we try to answer can be described as 
the problem of determining certain parameters (the age and 
mass of a star) by comparing models' predictions with em- 
pirical evidence (effective temperatures, luminosities, radii, 
dynamical masses of stars in binary systems). In order to 
do so we used a Bayesian approach, which allows to fully 
exploit the data. 

One of the main advantages of the Bayesian approach 
over the frequentist one is the possibility of using the avail- 
able information about the model parameters -the prior 
probability. Thanks to Bayes' Theorem this information is 
naturally included in the calculation of the new parameters' 
probability after additional evidence is collected -the poste- 
rior probability. In this way it is possible to further constrain 
the models' parameter space in an iterative process of re- 
finement. On the other hand the main disadvantage is that 
often the whole space of possible models is not accessible. 
This means that the normalizing factor appearing on the 
r.h.s of Bayes' Theorem -eq. ([l]) below-, can not always be 
evaluated. In such cases it is impossible to rigorously com- 
pute the normalized probability for a model to be correct 
given the empirical evidence. Nevertheless it is still possi- 
ble to compare and choose between two different models, by 
taking the ratio of the posterior probabilities thus removing 
the normalization factor. 

The notation we will adopt is the same as in 
|j0rgensen fc LindegrenI |200^ -hereafter JL05. The method 
described in JL05 has been successfully used to provide stel- 
lar ages and masses for the stars in the Geneva-Copenhagen 
survey (see iNordstrom et al.l |2004| ; iHolmberg etal] |2007| . 
iooi). It is a general method that can be applied to many 
other astrophysical cases, such ours. We retain most of the 
formalism of JL05, even though we customized the method 
for use with PMS objects. One major difference between the 
present work and JL05 is the use of the full covariance matrix 
when dealing with correlated variables such as luminosities 
and temperatures as determined for stars in eclipsing binary 
systems (see Sect. 12. 2| . An other important difference is the 
adoption of prior distributions that are appropriate for our 
observed sample. The flexibility of the Bayesian approach 
introduced by JL05 lies also in the opportunity of choosing 
the prior distributions that are best suited for the particu- 
lar problem. For example our objects have well-determined 
masses, hence we adopt these values as priors (see Sect. 12. 4] ). 

Let q be a set of observational quantities (or any combi- 
nation of them), for example temperature and luminosity or 
gravity and temperature. Let p be a set of model parameters 
and H a set of meta-parameters identifying a class of models. 
We introduce this distinction between p and H for practical 
reasons. The parameters p are the triple {t,ij,,(), i-e. age, 
mass and metallicity of the model. The meta-parameters 
H are instead the mixing-length parameter, a, the primor- 
dial helium abundance, Vp, and the helium-to- metals en- 
richment ratio AY/AZ. These three meta-parameters are 
chosen on the basis of some considerations and can be re- 
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garded as fixed inputs for the evolutionary models library 
as a whole, which gives them a different status compared to 
the p set. Th e a parameter is usu ally calibrated on a so- 
lar model fsee lBasu fc AntTall2008l . and references therein). 
The Fp value is constrained by big-bang nu cleosynthesis 
and observation of metal poor Hii regions (llzotov et al. 
20071: iPeimbert et"al] l2007l : IPunkley et aLll2009l : ISteigman 



2010r ). Ay/A.Z is constrained by chem i cal evolution mod- 



els o f the Galaxy (|Romano et al.l l2005l : ICarigi fc PeimberO 
I2OO8I ) or by comparing the absolute ma gnitude of unevolved 
nearby dwarf stars with stellar models ([Jimenez et al.|[20o3 : 
ICasagrande et al.ll2007l : iGennaro et al.ll2010l ). 

The H triple is usually fixed for any set of stellar tracks 
or isochrones available in the literature. Nevertheless there 
is no strong reason to assume that the solar calibration of 
a has to be suitable also for PMS stars of any mass (see 
e.g. the discussion in Montalban et al.ll2004l : iTognell i et al.l 
l201lD . Also Fp and AY/AZ are known with some uncer- 
tainty, which is quite large specially for the latter. Hence, 
having the opportunity to calculate our own stellar models, 
we allowed for variations of the meta-parameters, calculat- 
ing stellar models libraries for a total of nine combinations 
of them (see Sect. [3] for a detailed description). In principle 
other meta-parameters exist, such as the relative distribu- 
tion of metals (the mixture), the opacity tables, the equa- 
tion of state. Indeed any choice of physical inputs identifies a 
class of models, but we will not explore the effects of changes 
in the microphysics. 

In the following we will drop the H term and, if not 
explicitly stated, we will refer only to one particular class 
of models, i.e. one fixed choice of Hj = (oj, Ypj, AY/AZ^.). 
We will come back to the comparison of different classes in 
Sects, [m and El 



2.1 Definition 

Bayes' Theorem states that the posterior probability of the 
parameters p, given the observations q, is; 

/(«7|P)/(P) 



/(p|q) = 



/(q) 



(1) 



The probability f{q\p) of observing the quantities q given 
the parameters p is proportional to C(p\q), the Likelihood 
of the parameters p given the evidence q. The quantity f(p) 
is the prior distribution of the parameters, which incorpo- 
rates the information already available about them. The nor- 
malizing factor, /(q) = / f{q\p)f{p)Ap is called marginal 
distribution; it represents the probability of observing new 
evidence q under a complete set of mutually exclusive hy- 
pothesis, i.e. under all possible values for p. 

To calculate the integral, we should have access to all 
possible models, for all possible sets of parameters (and 
meta-parameters), or at least the subset of all plausible mod- 
els, i.e. models for which /(p) is not negligible. Even though 
the integration is not possible here, this is not a problem. 
As long as we are interested only in comparing different 
classes of models or estimating the most probable set of pa- 
rameters within a single class, this can be accomplished by 
taking probability ratios, hence removing the normalization. 
Having considered that, we can then rewrite the posterior 
probability as: 



2.2 The Likelihood function and the 2- variables 
covariance matrix 

Equation © is identical to eq. (3) in JL05, but in our work 
we extend the definition of Likelihood to the case of pairs 
of observables with non-zero covariance. This is particularly 
important -and often neglected- when the observables used 
to determine the stellar parameters are luminosity and tem- 
perature of stars in eclipsing binaries. Because of the way 
the two quantities are de rived, they are strongly correlated 
(see lMathieu et al] |2007)'). Let the vector of observables be 
a 2D vector: q — {x,y); the definition of Likelihood in the 
general case is: 



^(p|q) = 



1 



2n o-y \J\ — p 



X exp 



2(l-p2 



. [y{p)~yf 



2p [x{p) - x] [y{p) - y] 



(3) 



Here x and are the measured value for the observable 
X and its uncertainty, respectively (the same for y and ay). 

The quantity p — ^-^ — - is the correlation coefficient of x 

O'xO'y 

and y. The quantities x{p) and y{p) are the values predicted 
by the model for the parameter values p. 



2.3 The value of the covariance between 
luminosity and temperature 

In the case of eclipsing binaries the quantities that can be 
determined from the light curve are the effective tempera- 
tures ratio between the primary and secondary, T'eff,i/Teff,2, 
and the radii, Ri and R2. The temperature of the pri- 
mary has to be inferred by other indicators such as some 
temperature-sensitive lines in the spectrum and a subse- 
quent spectral type - temperature conversion. Luminosities 
are not directly measured, but derived using Stephan's Law: 
L = 47rasBi?^reff''. 

On the other hand, since the most used tool of stellar 
evolution is the Hertzsprung-Russel diagram (HR diagram) , 
where logL and logTcff are displayed, most of the analy- 
sis of binary systems is done in the HR diagram. So it is 
useful to have a proper treatment of the covariance matrix 
for luminosity and temperature. Given Stephan's law, the 
standard deviation of log L is calculated as: 



o^ioR L — \ -, — — - 2—— + 4- 



In 10 V -R Tef 



(4) 



fip\q) cc C{p\q)f{p) . 



(2) 



where an and ar^^^ are the uncertainties on radius and tem- 
peratur^. The covariance between log L and log Teff is given 
by: 

Cov(log L, log Toff) = 2 Cov(log R, log T^g) + 4 Var(log T,s) 
^ 4Var(logTeff) (5) 

where Var(logroff) = (yfosT^ir 

In both eqs. Q and ^ we assume that temperatures 



Note that In is the base-e logarithm and log the base-10 
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and radii have vanishing covariance. This is not necessarily 
true for echpsing binaries, for which they are derived from 
the same hght curve using fitting algorithms. However we 
are forced to neglect the corresponding term in the total 
luminosity-temperature covariance since we do not have ac- 
cess to the covariance matrix for radii and temperatures. 
Nevertheless we expect this covariance to be small specially 
when radii and temperatures are derived by multiple fitting 
of light-curves obtained independently in several photomet- 
ric bands. 



2.4 Prior distributions 

We will make use of different types of prior distributions for 
the parameters. Here is a brief description of each of them: 

Mass: for most of the systems in our sample the dynam- 
ical mass is the observable that is known with the best pre- 
cision. For this reason we will often use a Gaussian mass 
prior /(p) (X exp{ — i [(/i — mob)/o-m„bl^}- Apart from bet- 
ter constraining the mass values, this particular kind of prior 
is very informative and helps constraining also the stellar 
ages (see e.g. Sect.|4ll. 

Metallicity: for most of the systems measurements of 
[Fe/H] are available, too. These can also be included as 
priors after converting [Fe/H]c,b into the corresponding Z^b 
value. The value of Zoh depends on [Fe/H]ob but also on 
hYI l ^Z, Yp and (Z/X)c^ (for details see equations (1) and 
(2) in lGennaro et al.ll201(t ). As a consequence we use differ- 
ent Zoh values for a given [Fe/H]ob when we compare actual 
data to different classes of models S ^ . 

If we assume that [Fe/H] errors are distributed as Gaussians, 
the corresponding error distribution in Z is calculated as: 



/pHor(^) = GWZ)) 



d0(Z) 
dZ 



(6) 



where we introduced (l>{Z) = [Fe/H] and [Fe/H] is regarded 
as a function of Z; Q is a. Gaussian function. In the general 
case, the derived /prior (^) can be asymmetric and, in par- 
ticular, very different from a Gaussian. Nevertheless, given 
the typical errors in our data set (i.e. (T([Fe/H]) ~ 0.1 dex), 
the departure from Gaussianity is very small and can be 
neglected as can be seen in Fig. [2 where the Gaussian dis- 
tribution in Z and that obtained as in Eq. ((6]) are shown. 
In both cases, we start from an observed [Fe/H] distributed 
as a Gaussian and with ^ = 0.1 dex and a = 0.1 dex. The 
mean value and va riance of the G ( Z) are determined using 
eq. (1) and (2) of iGennaro et all (|2010l ) and simple error 
propagation rules. 

Moreover we can not always be sure that the error distri- 
bution on [Fe/H] is itself Gaussian. [Fe/H] errors certainly 
have a random, Gaussian component that can, however, 
be smaller than the overall uncertainty due to the poorly 
constrained systematics. With this in mind, we opted for a 
Gaussian functional form for the prior in Z. 

Age: in the case of two stars in the same system, co- 
evality might be considered as an additional prior, assuming 
that these two stars formed at the same time. In this case, 
we simply impose that both stars are coeval by multiplying 
their age marginal distributions (see below) hence getting a 
system age distribution. 
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Figure 1. Comparison between a Gaussian prior in Z and the 
prior that is derived as in eq. The assumed starting [Fe/H] 
distribution is a Gaussian with /i = 0.1 dex and cr = 0.1 dex. 



2.5 Marginal distributions, best values, 
uncertainties and relative precision 

As in JL05, the integration of eq.([2ll with respect to all pa- 
rameters but one, pi, yields the marginal distribution for pi. 
From this distribution it is possible to determine the most 
probable value for pi and its confidence interval. The two pa- 
rameters we are interested to determine are the stellar age 
r, and mass /i. We will use the same symbol as in JL05 for 
the age marginal distribution, G(r), and analogously define 
the mass marginal distribution, H{fj,). By writing explicitly 
the triple of parameters, we define: 



G(r) 



'C(-r,At,C|q)/(T,At,C)d/i dC 



Hifi) = / £(r,^,Clq)/(r,M,C)drdC 



(7a) 
(7b) 



JL05 demonstrated that the mode of the marginal dis- 
tribution is a more robust indicator than the mean for es- 
timating stellar ages. This is particularly true for strongly 
asymmetric distributions or distributions showing multiple 
peaks. We also adopted the mode as the best value estimator 
but changed the definition of the uncertainty interval with 
respect to JL05. If A is the total area under the distribution 
curve, F{x), we define the confidence interval [xmin, Xmax]: 



F{x)dx^ / F{x) 



dx = 0.16A 



(8) 



where we assume that the variable x is defined in the interval 
[a;i,a;u]. In this way 16% of the total probability is rejected 
on each side of the confidence interval. This definition co- 
incides with that of a la interval in the case of a Gaussian 
distribution. We followed again JL05 in the definition of the 
relative precision, e: 



XmsLX / Xn 



(9) 



Using this definition it is possible to compare the quality of 
different age and mass determinations. The worst relative 
precision is attained when the marginal distribution is fiat. 
In this case, assuming again x £ [xi, Xu] we have: 



100 v-^" 



Xl) 



Xl + 
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We calculated models in the age interval [0.5, 100] Myr and 
mass interval [0.2,3.6] Mq, hence the worst relative preci- 
sions attainable are e(r) ~ 1.26 and e(/i) « 1.03. 



2.6 Comparison of different classes of models 

Comparison of two classes of models is possible by calculat- 
ing the Bayes Factor, i.e. the ratio of the evidences for the 
two classes. The evidence itself is defined as the integral of 
the likelihood marginalized over the model parameters prior 
distributions. Hence the Bayes Factor for the i-th and j-th 
class of models is: 



where the evidence for each class is defined as: 



(10) 



/(qis) 



/(<7,p|S)dp 

,f(qb, s)/(p|E;)dp 



(11) 



The Bayes Factor tells us nothing about the best values 
of the parameters p, but it can be used to estimate which 
class of models -which set of meta-parameters H- gives an 
overall best-fit to the data. Strong deviations of BFij from 
one indicate that one class of models is a significantly better 
choice than the others. 



3 THE SET OF MODELS 

In the present analysis we used the very recent PMS tracks 
from the Pisa databas^EI which contains a very fine grid of 
models for 19 metallicity values between Z — 0.0002 and 
Z = 0.03, three different initial helium abundances and three 
values of the mixing-length parameter for each metallicity. 
The models have been computed using an updated version 
of the FRANEC evolutionary code which take s into account 
the s tate-of-art of all the input physics (see iTognelli et al.l 
I2OIII . for a detailed description). Here we briefly summarize 
the main characteristics of the code that are relevant for the 
present work and deeply affect both the morphology and the 
position of the PMS tracks in the HR diagram. 

We adopted the equation of stat e (EOS) released in 
2006 by the OPAL group (see e.g., [Rogers fc Navforiovl 
12002! ) . the OPAL high-temperature rad iative opacity re- 
leased by the same group in 2005 (see e.g..'lglesias fc Ro gers! 
[1996 1 for logr[K] > 4.5, and the jerguson e t al. (2005) low- 
temperature radiative opacities for logr[K] 4.5. The ra- 
diative opacity tables, both for low and high temperatures, 
are co mputed assuming the s olar-scaled heavy-element mix- 
ture bv lAsplund et al.! (|2005l '). 

The outer boundary conditions, required to integrate 
the stellar structure equations, have bee n taken from the de- 
tailed atmosphere models computed bylBrott & Hauschildt' 
l|2005! ) for Toff ^ 10 000 K, and by CasteUi fc Kurucz (2003) 
for higher temperatures. 

Convection is treated according to the Mixing Length 



^ The database is available at the URL: 
|http: //astro .df .ualpl ■ lt/st6llar-models/| 



Theory jBohm-Vitensd " l958|). f ollowing the formalism de- 
scribed in ICox fc Giulin f 1968! ). We used the classical 
Schwarzschild criterion to evaluate the borders of convec- 
tively unstable regions. 

The hydrogen bu rning reaction rate s are from the 
NACRE compilation (|Angulo et al.! Il999h . with the ex- 
ception of the "N(p,7)^'^0 from the LUNA collaboration 
I II ^ 

(|Imbriani et al. 2005). The code explicitly follows the chem- 
ical evolution of the light elements (D, "^He, Li, Be and B) 
from the early phases at the beginning of the Hayashi track. 
The models are evolved starting from a completely formed 
and fully convective structure, neglecting accretion. 

We extracted from the Pisa PMS database tracks for 12 
metallicities, namely Z = 0.007, 0.008, 0.009, 0.01, 0.0125, 
0.015, 0.0175, 0.02, 0.0225, 0.025, 0.0275, 0.03. The purpose 
was to cover the full range of metallicities for the observed 
sample of stars. For models in this range of metallicities 
we adopted an initial deuterium abundance Xp = 2 ■ 10~^ , 
suitable for pop. I stars (see e.g., Vidal-Madiar et al. l ll998l : 
iLinskv et aUl2006l : ISteigman et al.!l2007l ). 

For each value of Z the initial helium abundance, Y, 
has been obtained by the linear relation, 



Y = Yp + Z 



AY 

~Kz 



(12) 



where Yp is the primordial helium abundance, and AY/ A.Z 
is the helium-to-metals enrichment ratio. For Yp we adopted 
both the recent WMAP estimation Yp = 0.2485 (see 
e.g., ICvburt et al] |2004|; [Ste igmanI l200 6|) and a lower 
value Yp = 0.230 (iLeaueux et al., .1979: Pagel fc SimonsonI 
1 1989! : lOlive et al.] 119911 ). In the first case we used both 
AY/ AZ = 2 as commonly adopted in l i teratu r e (see 
e.g., iPagel fc Portinari! Il998l: Ijimenez et al] l2003l : iFlvnnI 
|2004 ICasagrande et al.! l2007h and AY/AZ = 5 that is 
the extreme value suggested by recent analysis (see e.g., 
ICennaro et al.![2010l . and references therein), while for Yp = 
0.230 we fixed AY/AZ = 2. Hence, for each value of Z, we 
computed models with three initial helium abundances. 

The efficiency of superadiabatic convection is 
parametrized by the a parameter where the mixing 
length I is given by ^ = aHp and Hp is the pressure 
scale-height. Following the usual procedure of calibrating 
the mixing length efficiency using the Solar observables, 
we obtained a = 1.68 for our reference set of models. 
However there is no strong reason to adopt this value 
for stars in different evolutionary phases compared to the 
Sun. Recent analysis of PMS stars in binary systems (see 
e.g., ISimon et all [2OO0I : ISteffen et all I2OOI! : IStassun et all 
I2OO4!) and studies of lithium depletion in young clusters 
( Ventura et"al]ll998l : [b'Antona fc Montalbanll2003l ) suggest 
a sub-solar efficiency of the superadiabatic convection in 
low-mass PMS stars. Therefore we decided to adopt tracks 
for three a values, i.e. a = 1.2 (low efficiency), a = 1.68 
(our solar calibrated), a = 1.9 (high efficiency). 

Taking into account the three Y values and the three a 
values, we computed models for each metallicity using a total 
of nine combinations of the H triples of meta-parameters. 

The tracks have been computed for a very fine grid of 
masses, with a spacing of 0.05 Mq in the range M = 0.2 — 
1.0 Mq, of O.IMq for the range M = 1.0 - 2.0 Mq and 
of 0.2Mq for the range M = 2.0 - 3.0AfQ. The tracks have 
been further interpolated on a finer mass-grid with a spacing 



6 M. Gennaro, P. G. Praia Moroni and E. Tognelli 



of 0.01 Mq and in age with a spacing of 0.05 Myr in the full 
mass range. This was done in order to achieve a very high 
precision in the determination of both the mass and age for 
the observed stars. 



4 SYNTHETIC DATA SETS: TESTING THE 
METHOD 

In order to check the accuracy of our method we first tested 
it against simulated data. As already demonstrated by JL05, 
the precision of inferred ages and masses is related to the 
detailed morphology of isochrones and tracks in the HR di- 
agram. Depending on the mass and the age of a star, its 
evolution might be faster in some parts of the diagram than 
in others. The evolutionary speed of a star along its track 
together with the positional uncertainty in the HR diagram 
determine the absolute precision of the method. Although 
the PMS is globally a very fast evolutionary phase com- 
pared to the MS, there are some stages -i.e. the descent of 
the Hayashi tracks- slower than others -i.e. the approach 
towards MS along Ifeyney tracks when radiative cores are 
developed. Another general rule is that more massive stars 
evolve faster, leading to progressively larger spacing between 
isochrones, and consequently a better relative precision in 
age determinations for a given evolutionary phase, as the 
mass increases. 

The evolutionary speed is inversely proportional to the 
age-gradient calculated along one evolutionary track. In the 
regions of the HR diagram where the evolutionary speed is 
large, the age-gradient is small and the precision in estimat- 
ing stellar ages is good. To explain this let us define an effec- 
tive temperature-luminosity error box, as given by observa- 
tional uncertainties. If we move it across the HR diagram, in 
regions with small age-gradient we will encircle models with 
similar ages, leading to a precise age determination, while in 
regions of high age-gradient (low evolutionary speed) within 
the same box there will be models with very different ages, 
leading to a less precise age estimate. In Fig. [2] (left panel) 
we show the evolutionary speed. Moving from blue-purple 
towards red-orange regions the evolutionary speed increases 
allowing for progressively more precise age determinations. 

In analogy to the age-gradient along a track we can 
calculate the mass-gradient along an isochrone. In this case, 
regions of small mass-gradient are regions where the stellar 
masses can be estimated with better precision. The mass- 
gradient along isochrones is shown in Fig. [21 right panel. In 
this case regions with lower mass-gradient (better precision) 
are in red. 

To explore how the position in the HR diagram affects 
age and mass determinations, we generated a sample of syn- 
thetic eclipsing binaries (EBs) for different masses and ages 
(see Table [TJ. 

The simulations of EBs have been done by selecting 
models from the PISA stellar library with a — 1.68, Yp — 
0.2485 and AY/AZ = 2. We fixed the Z value to Z = 
0.0125, similar to our solar model {Zq = 0.0137). For 
each combination of masses and ages we generated 100 
systems. We added random Gaussian uncertainties to the 
quantities predicted by the models using standard devi- 
ations values equal to the typical errors in our data set 
(cTM = 0.015-0.020 Mq, criogTrff = 0.015 dex, an = 0.057?© 



Table 1. Mass and ages combination for the simulated EBs. 
Systems C.x have the same masses of the binary system RXJ 
0529.4-1-0041 A. 



System 


Primary Mass 


Secondary Mass 


Age 


[Mq] 


M[Mq] 


[Myr] 


A.2 


2.0 


0.5 


2 


A.5 


2.0 


0.5 


5 


A.8 


2.0 


0.5 


8 


B.2 


1.0 


0.3 


2 


B.5 


1.0 


0.3 


5 


B.8 


1.0 


0.3 


8 


C.5 


1.27 


0.93 


5 


C.IO 


1.27 


0.93 


10 


C.15 


1.27 


0.93 


15 



and (JiogL — 0.1 dex). To simulate the observed errors be- 
haviour of EBs, we allowed for random errors in the pri- 
mary star temperature. Toff,!, keeping the ratio of primary- 
to-secondary effective temperatures, T'eff,i/T'efT,2, fixed. We 
allowed for independent errors in the radii, Ri and J?2. Lu- 
minosities are calculated from the temperatures and radii 
after the errors have been added. 

For each simulated system, we applied our Bayesian 
method to recover the best ages and masses. Since we have 
chosen to fix Z = 0.0125, we fixed it also in the recov- 
ery method, which is equivalent to using a prior /(^) = 
5(C — 0.0125) in equations (ffa)l and (l7b)l . We ran the method 
both with a flat prior defined over the whole mass inter- 
val and applying a Gaussian prior on the simulated mass. 
The Gaussian prior is centred on the simulated value of the 
mass and its a is of the order of the typical error for the 
dynamical masses available in the literature (few %). We 
obtained stellar ages for the single stars and also for the 
systems. In this last case coevality is imposed by consider- 
ing Gc(t) — Gp{t) X Gs{t), i.e. the product of primary and 
secondary marginal age distributions. 

The results for the complete set of simulations are 
shown in Table [2] where we report the percentage of cases 
in which the simulated age and mass fall within the confi- 
dence interval. From the table it is clear that the method 
is very successful in recovering the simulated values. When 
no systematic errors are present (in both the model and the 
data) and if the random error estimates are reliable, we can 
expect the results of the method for real data to be very 
robust. Unfortunately this ideal situation is seldom realized 
in reality, but it is worth noticing that the method is intrin- 
sically able to give a good fit for almost all the regions of 
the HR diagram. 

The fraction of good recoveries or success rate can be 
related to the position of the stars in Fig. [2] For example 
the secondary star of the B.x systems -a 0.3 Mq star- is 
moving towards slower phases of its Hayashi tracks and con- 
sequently the fraction of good age recoveries is decreasing 
with increasing age. Several other things are worth noticing 
about the recovery fractions in Table [2l 

• First, we note that the actual number of good recoveries 
is a complex function of the stellar position in the HR dia- 
gram. It is true that the mass- and age- gradient visualization 
of Fig.[2]can help understanding this function. On the other 
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Figure 2. Left: evolutionary speed along stellar tracks in the HR diagram. The age-gradient is inversely proportional to this quantity. 
The colour coding is a scale from red-orange (fast evolution, small age-gradient, good age determination) to blue-purple (slow evolution, 
large age-gradient, bad age determination). Right: mass-gradient calculated along isochrones. Red-orange regions are regions of low mass- 
gradient (high precision in mass determination) , blue-purple regions are regions of high mass-gradient where masses are determined with 
worse precision. In both panels symbols indicate the positions of the simulated stars before the random errors are added. Superimposed 
in white are some reference tracks and isochrones. 



hand we have to warn the reader that what we indicated 
as mass- and age-gradients are partial derivatives calculated 
along isochrones and tracks respectively. Therefore, they do 
not fully represent the real gradient. As a consequence there 
are regions of the HR diagram where the recovery fraction 
is different from what one might naively expect by looking 
at Fig. [2] alone. For example, the mass recovery fraction of 
the primary star of the B.8 case is lower than the B.5 case, 
even though the B.8 case is in a zone of lower mass-gradient 
(better mass resolution). 

• Secondly, we want to mention the power of the coeval- 
ity prior. Fig. |3] shows on the left the distribution of the 
difference between the logarithm of the best fit age and the 
logarithm of the simulated age for each star. The standard 
deviation of this distribution is cr = 0.185 dex. The inset 
plot shows the recovered ages for those stars giving a bad 
fit, meaning that the simulated age is outside the 68% con- 
fidence interval. The fraction of these bad-fit cases is 19% of 
the total simulated stars. 

In the central panel the stars in each system are paired to- 
gether and the age of the system is evaluated from the com- 
posite age distribution Gcij). The differences between the 
resulting best ages and the simulated ages show a much nar- 
rower distribution, with a — 0.062, almost 3 times smaller 
than the a for the single stellar ages. Hence using the coeval- 
ity prior strongly reduces the error in the best age estimate. 

The rightmost panel of Fig.[3]shows the difference between 
the logarithms of the primary and secondary components' 
ages for each pair. This distribution has a, a — 0.257. We 
note that in 13% of the cases the two components are found 
to be non coeval, meaning that the two uncertainty intervals 
do not overlap; these systems are shown in the inset diagram. 
The fraction of non coeval recoveries varies strongly among 
the simulated systems and the results are summarized in 
Table O The simulated binary that causes most of the non 
coeval fake detections is the system A. 8. Out of 100 simu- 



lated cases, 95 give non coeval results. This is due to the 
fact that the primary star of the 100 simulated A. 8 systems 
is very close to the main sequence. Therefore G(r) for such 
a star is a very flat function, with a broad uncertainty in- 
terval. Moreover the position of the peak of the distribution 
strongly depends on where the simulated system is scattered 
once that the observational errors are added. On the other 
hand, the age of the secondary is well recovered, with a nar- 
row uncertainty interval. As a consequence, in the 95% of 
the cases, the two intervals become disjoint and the two stars 
are considered as non coeval. 

The fraction of non coeval recoveries is a complex func- 
tion of the position in the HR diagram of both stars. For an 
observed sample of binaries, the fraction of total fake non 
coeval systems to be expected depends on the actual distri- 
bution of relative positions of primaries and secondaries in 
the HR diagram. 

Nevertheless we point out that the tail of non coeval sys- 
tems disappears in the central diagram of Fig. [3] i.e. no 
system age is found to be in disagreement with the sim- 
ulated age, even when the two components were formally 
non coeval. In our simulated cases, the non coevality is 
mostly due to one of the stars being in a region of very large 
age-gradient, e.g. close to the main sequence. On the other 
hand the companions are generally in a region of smaller 
gradient. Therefore, while the age of one star is badly de- 
termined, the other component still has a very informative 
G(r), which drives the composite distribution towards the 
simulated value. This fact can be used to determine the age 
of main sequence stars in systems where the companion is 
still on the PMS, for example. 

In addition, for the full sample of simulated cases, the use 
of the composite age distribution allows to find a system age 
that is in better agreement with the simulated one. The a of 
the distribution of recovered minus simulated age is reduced 
when the coevality prior is imposed. 
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Table 2. Percentage of cases in which simulated ages and masses 
are recovered within the 68% uncertainty interval. P indicates 
the primary stars in the systems, S the secondary, C stands for 
Coeval, indicating the cases in which the product of the marginal 
age distributions of the primary and secondary star Gp (t) X Gg (r) 
is used to infer the age of the whole system. Flat and Gaussian 
are the adopted mass priors. 



System 




A.2 






A. 5 






A. 8 




Flat 


P 


S 


C 


P 


S 


C 


P 


S 


c 
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83 
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74 
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Masses 
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68 




98 
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78 


79 
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B.5 
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94 
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96 
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100 
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System 




C.5 






C.IO 






C.15 




Flat 


p 


S 


C 


P 


S 


C 


P 


S 


C 


Ages 


95 


98 


89 


83 


93 


78 


87 


86 


94 


Masses 


82 


74 
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68 




86 


84 
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94 


93 


72 


88 


58 


71 


77 


78 
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100 


100 




100 


98 




100 


100 





These results are very important. It is sometimes noted 
in the literature that evolutionary models are not able to fit 
binary data for the same age and ad hoc solutions are in- 
voked to reconcile the models and the observation. Here, we 
demonstrated that not being able to reproduce both com- 
ponents in a binary system with the same isochrone does 
not necessarily imply that the stars are not coeval or that 
they are coeval but the models are not able to reproduce 
this coevality. On the contrary, an age difference or even an 
age mismatch can simply be a consequence of the observa- 
tional errors in the HR diagram. It is indeed possible that 
the random scatter of the positions of coeval stars acts in 
opposite directions for the two components, making one look 
older and the other younger, to the point that they might be 
considered as non coeval. The actual expected artificial age 
difference depends on the region of the HR diagram where 
the two stars are located, and on their errors. 

As an example we consider the work by 
iKraus fc Hillenbrand! (|2009l ) who analyse the binary 
population in the Taurus- Auriga association. They found 
that, in general, stars in physical pairs are more coeval 



Table 3. Fraction of non coeval fake recoveries for the simulated 
systems. 



System 


Simulations 


Non coeval recoveries 


A.2 


100 


% 


A.5 


100 


10 % 


A.8 


100 


95 % 


B.2 


100 


% 


B.5 


100 


% 


B.8 


100 


8 % 


C.5 


100 


% 


C.IO 


100 


% 


C.15 


100 


2 % 


Total 


900 


13 % 



than the association as a whole, with significantly smaller 
intra-binary age spread than for randomly paired stars 
selected among the association's members. Nevertheless 
the authors also find that some of the binaries show an 
intrinsic age spread larger than what observed for the bulk 
of the pairs and not consistent -within the errors- with the 
hypothesis of coeval pairs. They suggest that these outliers 
can be multiple system with unrecognized companions 
or stars seen in scattered light or also stars with disk 
contamination. While this can certainly be the case, we 
want to point out that the observational errors themselves 
-even when one might think that they are completely under 
control as in our simulations- can be partially responsible 
for an artificially large age-spread (or non coevality) within 
a binary system. The fraction of objects for which this fake 
non coevality might be expected is a complex function of 
the actual distribution of observed stars in the HR diagram, 
i.e. of their mass, mass-ratios and ages. 

• As a third remark we want to emphasize the power of 
the Gaussian mass prior, i.e. its large informative value for 
the recovery of stellar masses. In almost all cases where the 
Gaussian mass prior is imposed, the recovery fraction raises 
to 100%. This might look obvious but we have to recall that 
dynamical masses are usually the most reliable data avail- 
able for a binary system. The impact of the Gaussian prior 
is less strong when the results for stellar ages are compared. 
Among the 27 cases displayed in Table [21 14 favor the flat 
prior, while in 12 cases the Gaussian prior gives better re- 
sults. The primary star of the A.8 system gives equal frac- 
tions of 0% recoveries for Gaussian and flat prior (the star 
is on the MS). 

We investigated this behaviour and observed that the im- 
pact of the Gaussian prior actually depends on the detailed 
morphology of the 2D posterior probability in the — t 
space. As an example, consider the primary star for the C.IO 
case. The success-ratio in the age recovery decreases from 
83% to 72% for this star when the Gaussian mass prior is 
applied. 

We show in Fig.Uone of these drop-out cases. The reason 
why the age for this particular star is not recovered any- 
more when the prior is included can be easily understood. 
The Gaussian prior is causing an increase in the relative 
precision for the marginal distribution, GGauss(T). This can 
be seen in the right panel of the figure where the GGauss(T) 
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Figure 3. Left: distribution of the difference between the logarithm of the recovered age for the single stars and the logarithm of 
the simulated age. The inset shows the distribution of age differences only for the stars for which the simulated age is outside of the 
68% confidence interval Centre: same as left, but the recovered age is obtained from the composite age distribution for each simulated 
binary system. Right: Distribution of the differences between the logarithm of the primary's recovered age and the logarithm of the 
secondary's recovered age for each system. The inset shows the distribution of age differences for those systems whose primary and 
secondary components have disjoint 68% age confidence interval and are considered as non coeval. 



distribution has a clearly larger mode and is narrower than 
Gmaxir). One can notice that the right border of the con- 
fidence interval for G'Fiat(T) is already quite close to the 
simulated age value of 10 Myr, which is barely within the 
68% confidence interval. Hence it is the shrinking of the con- 
fidence interval -when the Gaussian prior is applied- that 
causes the simulated age value to drop out of the 68% confi- 
dence interval. Also the actual number of dropouts is related 
to the detailed structure of the 2D mass- and age-gradients. 



5 THE DATA SET 

The number of PMS stars with direct mass measurements 
amounts nowadays to about 30 objects. The sample we used 
consists of 25 PMS and 2 MS stars, whose properties are 
summarized in Table [l] 

Among the 27 objects 10 are PMS stars found in eclips- 
ing binary systems (EB). We included the MS 3.1 A/q star 
TYCrA A and 2.1 Mq star EK Cep A in our sample as well, 
since their companions are PMS stars in EB systems and we 
tried to fit both components of binary systems when possible 
(see Sect.[8|. For EBs the masses, radii, and effective tem- 
perature ratios of the two components can be accurately in- 
ferred. However the absolute values of the effective tempera- 
tures, that rely on the determination of the primary's spec- 
tral type and some spectral type-Tefi relationship, can be 
affected by systematic errors a nd constitute a severe source 
of uncertainty. As an example iHillenbrand k, White! (|2004l ) 
pointed out that the main sequence empirical scales used for 
deriving temperatures from the spectral types might cause 
systematic temperature offsets when applied to PMS stars, 
due to the differen t values of surface gra vity for a given spec- 
tral type (see also lLuhman et al]|l997l 'l. 

Of the 27 stars in the sample, 6 are found in astromet- 
ric/spectroscopic systems (AS), i.e. systems in which the 



components can be resolved as separate point sources us- 
ing interferometry. Combining astrometry and line-of-sight 
velocity measurements the masses of the components are de- 
termined in a distance- independent way. With this technique 
the radii of the components are not measurable, though. 

The last 9 objects have masses measured using their 
circumstellar disks keplerian velocities obtained by means 
of spectroscopy. The mass for the central star can be deter- 
mined only if the linear value of the orbital radius at which 
the velocity is measured is known. Hence stellar masses are 
in this case distance-dependent. The two stars in the UZ 
Tau E system form a binary, and their masses are separated 
using combined spectroscopic measurements for the circum- 
stellar disk and the stellar velocities (DKS). 

A very s imila r sa mple was alr eady studied in 
iMathieu et akl (|2007l ') and IStassunI (|2008l ') and we refer the 
reader to the first of these papers for a detailed description of 
the different observational techniques and the different kind 
of unc ertainties affecting them. Compared to iMathieu et al.l 
(20071) there are some distinct o bjects in our sample though. 
The 2M0535-5 brown dwarf EB l|Stassun et al.ll2006l 'l was ex- 
cluded because the stellar dynamical masses are smaller than 
those currently present in the Pisa databas e. The recently 
disco vered PMS EB ASAS J052821+0338.5 (|Stempels et al.l 
l2008h has been added t o the sample. We al so included the 
AS binary HD 113449 l|Cusano et all I2OIOI ') for which we 
have slightly different parameters from an updated analysis 
(Cusano, private communication). 

The luminosities and effective temperatures for our 
complete data set are displayed in Fig. [S] Overplotted are 
stellar tracks for a = 1.68, AY/AZ = 2 and Fp = 0.2485. 
We will refer to this set of parameters H as our standard or 
reference set. In the case of Fig. [S] the tracks are calculated 
for a value of Z — 0.0125, similar to the metallicity of our 
standard solar model (i.e., Zq = 0.0137). 
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Figure 4. A particular realization for the primary star of the C.IO system. In this particular case the age is well recovered in the flat mass 
prior case and not recovered when the Gaussian mass prior is imposed. Left: the HR diagram position for the star with superimposed 
some reference isochrones (9,10,11 Myr) and tracks (1.17,1.27 and 1.37 Mq). The cross indicates the position after errors are added, the 
filled circle represents the original position of a 1.27M0, 10 Myr star. Centre: 2D posterior probability contours. Dotted line is for the 
case with flat mass prior and solid line for the case with Gaussian mass prior. Right: Marginal age distributions. Vertical lines indicate 
the confidence interval. Dotted and solid line as for the central panel. 



For several systems, [Fe/H] values are available from 
direct spectrosc opic measurern e nts. W e used the [Fe/H] de- 
termin ations bv lD'Orazi et al] ([2009') and bv lD'Orazi et al] 
l|201lll for the systems in Orion and Taurus- Auriga, respec- 
tively. Only 4 stars are left without a [Fe/H] measurement. 
In order to convert the observed [ Fe/H] into the global 
metallicity Z^h, we followed eq. ( 2) in Gennaro et al. _ l|201(l ). 
adopting {Z/X)q = 0.0181 bv lAsplund' et al ] (|2009l ). Al- 
tho ugh the Pisa PMS mo dels have been computed adopting 
the lAsplund et ahl l|2005f ) heavy elements solar mixture, the 
seeming in consistency is in consequential since models com- 
puted with lAsplund et al. I i2005) and Asplund ct al. (2009) 
but with the same total metallicity Z are essentially indis - 
tinguishable (see detailed discussion in iTognelli et al.|[2"oilf ). 

A point worth mentioning is that among the objects 
listed in Table 3] there are a few cases that appear peculiar. 
Their location in the HR diagram is indeed incompatible 
with that of stars of similar mass. By looking at the available 
dynamical masses, luminosities and temperatures we iden- 
tified four of these peculiar objects, namely NTT 045251 A, 
UZ Tau Ea, BP Tau, and MWC 480 (see Table H). 

• NTT 045251 A is close to RXJ 0529.4 Ab and V1174 
Ori A in the HR diagram, with approximately the same 
luminosity and T^s. However NTT 045251 A is more massive 
than the other two stars by 0.4-0.5 Mq. We checked that 
this discrepancy can not be reconciled even by assuming that 
the metallicity of NTT 045251 A is 0.4 dex larger than that 
of RXJ 0529.4 Ab or VI 174 Ori A. 

• UZ Tau Ea, M = 1.O16M0, has a Tch similar to V1174 
Ori B, DM Tau, CY Tau, and NTT 045251 B which have 
lower masses, between 0.55-0.8 Mq. Moreover the star is 
colder than RXJ 0529.4 Ab and V1174 Ori A by about 850 
K in spite of their similar masses, luminosities, and [Fe/H]. 

• BP Tau, M = 1.320Mq, is significantly colder (ATeff > 
1000 K) and fainter (AlogL/L© ^ 0.9) than RXJ 0529.4 



Aa and EK Cep B, although it is slightly more massive and 
the metallicities are similar. In addition this star is fainter 
and colder than the 0.96 Mq HD 113449, which has a similar 
metallicity. This is not easy to explain because the minimum 
luminosity of a 1.3 Mq model is always larger than that 
achieved by a 0.9-1.0 Mq star approaching the ZAMS, as in 
the case of HD 113449 A. 

• MWC 480 (M = 1.65Mq and [Fe/H] = -0.01) is located 
in the HR diagram between EK Cep A (M = 2.02Mq and 
[Fe/H] = 0.07) and the RS Cha system (M = 1.87, 1.89Mq 
and [Fe/H] — 0.17). We checked with our models that the 
difference in [Fe/H] can justify neither the similar luminosi- 
ties of MWC 480 and EK Cep nor the higher luminosity of 
MCW 480 with respect to the RS Cha stars. 



6 THEORETICAL VS. DYNAMICAL MASSES 
- THE STANDARD SET OF MODELS 

In this section we show the comparison of the whole data 
set with our evolutionary models. We limited the analysis 
only to the standard set of models, i.e the H class with a = 
1.68,Ay/AZ = 2 and Yp = 0.2485. Section □ is dedicated 
to the comparison of models with different meta-parameters, 
H. Since radii are not available for each star in the data set, 
the comparison was done in the HR diagram. With the no- 
tation of Sect.[2l this means q — (logTcfi, log L/ Lq). In the 
case of stars with available [Fe/H] measurements, a Gaus- 
sian prior on the metallicity was applied after converting 
the [Fe/H] values and their errors into Z values with corre- 
sponding errors az'- 



/(C) = 



X cxp 



2a| 



In the other cases a flat prior for Z G [0.007, 0.03] was used. 
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Table 4. List of stellar properties. 
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The outcomes of the full data set comparison are sum- 
marized in Fig. [S] The stars have been divided into sub- 
groups: eclipsing binaries (EB), astrometric/spectroscopic 
binaries (AS) or disk kinematic stars (DK). The stars in 
the UZ Tau E system (DKS, i.e. disk kinematics plus spec- 
troscopy to disentangle the components) are included in the 
DK sample for simplicity. The three subgroups are displayed 
from left to right. Each panel shows a comparison of the rel- 
ative difference between model-inferred mass -Mmod- and 
measured dynamical mass -Mdyn . In the upper panels Mmod 
is derived by applying the Gaussian metallicity prior, when 
available. The lower panels show a comparison between the 
Gaussian Z prior case (empty symbols) and the flat Z prior 
case (full symbols). The symbols indicate the mode of the 
posterior probability; the asymmetric error bars indicate the 
68% confidence interval as described in Sect. (2] In the fig- 
ure the dynamical mass errors are not added to the error 
budget. This is meant to purely show the precision of the 
masses estimated from the models given the observational 



uncertainties. However this is not a bad approximation for 
the total mass error budget, given that the quoted errors are 
of the order of 1% for most of the dynamical masses and up 
to 10% only in very few cases (see Table [Jjl. 

By looking at the figure, it is clear that EB masses are 
well recovered in almost every case but for 'V1174 Ori B. We 
will discuss this particular object more in detail in Sect. 18.31 
The general agreement becomes progressively worse for AS 
binaries and DK stars. The AS binaries have masses that 
are underestimated by ~ 20 — 30% on average. The worst 
case is NTT 045251+3016 B, whose mass is underestimated 
by ~ 50% (see Sect. [8^ . For DK stars the trend becomes 
more negative with underestimates as low as 70% for UZ 
Tau Ea. Regarding the DK group, we have to point out that 
for the latter class of objects the uncertainty on the distance 
is not included in the dynamical mass error estimate. This 
uncertainty propagates linearly in the mass uncertainty and 
quadratically in the luminosity uncertainty. The DK stars in 
our sample stars are part of the Taurus- Auriga star-forming 
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Figure 5. The HR diagram for our data set. Tiie labels cor- 
respond to the ID column in Table |4] Superimposed are stellar 
tracks calculated with a = 1.68, AY/ AZ = 2 and Yp = 0.2485 
and Z = 0.0125. The values of the mass -in solar units- are dis- 
played on the left of the corresponding track. 

complex which is located at about 150 pc from the S un and 
has a radius of about 15 pc (see e.g lTorres et al.|[2009l ). Using 
an average distance to each star instead of its real distance 
may then cause a systematic error on the mass estimate of 
about 10% and on the luminosity of up to 20%. Part of 
the disagreement might also arise from the fact that the DK 
objects are T-Tauri stars, intrinsically variable. The temper- 
atures and luminosities we a dopt for them are all derived by 
iHillenbrand fc White! (|2004 ). The authors try to mmimize 
the effects of accretion luminosity using the Ic band to esti- 
mate the stellar luminosity. Nevertheless the estimated stel- 
lar luminosities might still be offset from their real values. 
An other problem might affect the temperature determina- 
tion. Temperatures are determined fr om spectral types us- 
ing re lations calibrated on dwarfs. As IHillenbrand fc Whitel 
l|2004l point out, there might be a systematic temperature 
underestimate due to the fact that PMS stars of a given 
spectral type are generally warmer than dwarf counterparts. 
If the real stellar temperatures would be higher part of the 
discrepancy in our mass estimates would be removed, since 
larger masses would be needed to reproduce the observed 
stellar properties. 

The lower panels of fig.[6]show a quite surprising result. 
In the cases where we applied the Gaussian prior on the 
metallicity the final results for Mmod are usually in worse 
agreement with the Mdyn values when compared to the flat- 
prior case. It appears that applying the metallicity Gaussian 
priors we obtain, in general, lower masses values. This be- 
haviour suggests that the Zoh values used here might be 
too low. If we force the metallicity to assume systematically 
lower values than the real values by means of the Gaussian 
prior, then we naturally obtain lower values for the best 
masses. This is due to the fact that the lower the metallicity 
used in a stellar model, the hotter and more luminous the 
model is for a given mass. 

We speculate here that part of the problem with the 
metallicity prior might reside in the low value of {Z/X)q 
used in the present work to convert the observed [Fe/H] into 



Z. In the recent years this value has undergone a drastic 
change, specially after the introduction of non-LTE and 3D 
hydrodynamical atmospheric models for the analysis of Solar 
abundances (see iAsplund et al.li2009l . for a review about the 
topic). The traditional value of Zq ~ 0.02 has been strongly 
revised t owards much lower v alues, down to Zq 0.013-^0.014 
(see e.g.. lSerenelli et al]|2009l . and references therein). There 
is still ongoing debate regarding solar heavy elements rela- 
tive abundances and total metallicity, and the uncertainty 
on the absolute values is still large. Nevertheless the most re- 
cent results are going in the increasing direction for {Z/X)q, 
thus reducing the difference with the traditional estimates. 
A change in {Z/X)q will be reflected directly into Zoh since 
(Z/X)oh = {Z/X)q X lo''^"/^' and an increase in this quan- 
tity would naturally lead to a systematically larger Mmod- 



7 ANALYSIS OF THE DATA USING 
DIFFERENT CLASSES OF MODELS 

As described in Sect. 12.61 the ratio of the evidences for 
two classes of models -the Bayes Factor, BF- can be used 
to quantify which class of model is better in reproducing 
the data. Significantly better evidence of a model over an 
other is claimed when BF < 0.1 or BF > 10, i.e. when 
the two evidences diffe r by one order of magnitude or more 
(IKass fc Raft,ervl[l995l ). Since the evidence is calculated by 
marginalizing the posterior distribution over all the param- 
eters of the model, the prior distribution of the parameters 
has to be consider ed as part of the model as well (see e.g. 
iBailer- Jonesll201ll. for an application of the Bayes Factor 
to discriminate between distinct models in a different astro- 
physical context). 

We calculated the evidence for each star in the 9 differ- 
ent meta-parameters cases and for each set H we considered 
4 different combinations of the prior distributions for the 
masses and the metallicities. The four combinations are: 1) 
flat mass prior - flat metallicity prior, 2) Gaussian mass prior 
- flat metallicity prior, 3) flat mass prior - Gaussian metal- 
licity prior, 4) Gaussian mass prior - Gaussian metallicity 
prior. Hence we actually have 9 x 4 = 36 classes of models. 
After evidences have been calculated in the 36 cases for each 
of the 27 stars of the sample, the BFs have all been calcu- 
lated by dividing each evidence value for a given star by the 
evidence value for that same star obtained using our stan- 
dard class of models. The latter is identifled by Yp = 0.2485, 
AY/AZ = 2 and Q = 1.68 for the 1) case. 

The numerical values of the BFs are reported in Ap- 
pendix S The Tables [M\ IX2l [Ml and [HI displav the BF 
values for the 1), 2), 3), and 4) prior cases respectively. Look- 
ing at each table, it is possible to see the change in the evi- 
dence among the different Hj sets, within one of the 4 prior 
cases. From one table to the other, the corresponding entries 
are calculated for the same H set but in the 4 different prior 
cases; by looking at these corresponding entries it is possible 
to see the role played by the prior choice and, actually, un- 
derstand whether this choice is leading to an improvement 
of the overall flt or not. 

The best way to compare different classes using the 
whole data set is to calculate the composite evidence for 
the full set of data. The natural extension of equation (1111) 
is that the evidence for the whole data set -represented here 
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Figure 6. Inferred masses from our standard set of isochrones compared with tiie dynamical masses from literature. Numbers correspond 
to the ID column of Table |4] From left to right : stars are divided in subset of eclipsing binaries (EB), astrometric/spectroscopic binaries 
(AS) and stars with masses inferred from disk kinematics (DK). Upper panels: best values are indicated by empty symbols. The asymmetric 
error bars represent the 68% confidence interval, as defined in Sect. [2] In these cases a Gaussian metallicity prior was applied when [Fe/H] 
measurements were available. Lower panels: Comparison of the best masses inferred when the metallicity Gaussian prior is imposed (empty 
symbols, same as upper panels) and when we use a flat prior even for stars with [Fe/H] measurements (filled symbols). 



by the set of observables {q}~ can be written as: 

k •' 

= n / /(q.|P,3)/fc(p|H)dp (13) 

k ■' 

The index k runs over the stars in the sample and fk in- 
dicates the specific prior distribution applied for the fc-th 
star. 

We multiplied the evidences of the 27 stars of the sam- 
ple to understand which class of model gives the best gen- 
eral result. In addition we restricted the product to stars 
belonging only to one type of system (EB, AS, DK). When 
using the whole data set we can compare only the 1) and 
2) cases because the Gaussian metallicity prior can not be 
applied to the 4 stars for which [Fe/H] measurements are 
not available. To compare also the 3) and 4) cases we addi- 
tionally restricted the sub-samples to only stars with [Fe/H] 
estimates. 

In almost all the cases, the class of models with 



Ay/AZ = 2,Yp = 0.23 and a = 1.2 meta-parameters is the 
one with the largest evidence. The only exception is the case 
of AS subset when only stars with known [Fe/H] are con- 
sidered. For this subset the largest evidence is attained for 
the AF/AZ = 2, Yp = 0.23 and a = 1.9 meta-parameters 
values. We have to point out that in this particular case only 
4 stars are part of the subset and the evidence is only 1.13 
times larger than in the AF/AZ = 2, Yp = 0.23, a = 1.2 
case. 

These results are the counterpart of what we have al- 
ready observed in Sect. [6] regarding the mass underestima- 
tion by standard models. It was clear -specially by looking 
at the upper panels of Fig. [6]- that the general trend for the 
standard set of models is to predict too low masses com- 
pared to Mdyn . A sim ilar trend can be observed in Fig. 3 of 
iMathieu et all l|2007D : here all the considered sets of mod- 
els show the same kind of behaviour in predicting too low 
stellar masses, with a mean difference that can be of the or- 
der of 20% or more. This suggests that the standard tracks 
are too hot and luminous when displayed in the HR diagram 
compared to the observed temperatures and luminosities for 
the given dynamical masses. A natural way to get a better 
agreement with observations is to use colder and fainter set 
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of models by adopting both a lower helium initial abundance 
and mixing length parameter a value. 

This explains why the best overall evidence is achieved 
by the set with a = 1.2, Yp = 0.230, and AY/AZ = 2. 
Similar low helium content for a given metallicity Z could 
be obtained also adopting the currently accepted primor- 
dial helium value, Yp = 0.2485, together with a very small 
helium-to-metals enrichment ratio {/SY/AZ < 1, for Z ~ 
0.01 0.02). However, both choices are quite unlikely. In 
the former case the Yp = 0.23 value is significantly lower 
than the recent independent results from extr agalactic Hll 
regions and Big Bang Nucl e osynthesis theorv (p^zotov et al 



regions and Jaig Jaang IN ucl e osynthesis theorv (lizotov et al. 
2007 : iPeimbert et al T l2007l : iDunklev et a' . I l2009l : ISteigmanl 



20101). Regarding the latter case, a Ay/AZ < 1 is smaller 



than the val ue suggested by both Galactic chemical evolu - 
tion models (|Romano et al.ll2005l : [Car igi fc Peimbert"2008'') 
and by nearby dwarf stars analysis (Jimene z ct al. 2003; 
ICasagrande et al.ll2007l : iGennaro et al.ll2010f) . 

We point out that when a single star is considered the 
value of the evidence for different classes of models are 
mostly of the same order of magnitude (see the Tables in 
Appendix 13 . Hence these global results are more sensitive 
to the few objects for which we observe major changes in 
the evidence between different classes. Nevertheless a gen- 
eral analysis is still important to understand the overall be- 
haviour of stellar models in comparison to available data. 

We also report that, for the whole sample, as for any 
subset of stars (EB, AS or DK), the largest value of the ev- 
idence is reached in cases 2) or 4), i.e. when the Gaussian 
mass prior is imposed. This is expected because the prior on 
the dynamical mass is a much more informative prior than 
a flat one defined across the entire mass range of simulated 
models (i.e. 0.2-3.6 Mq). We already pointed out in Sect. [3] 
that imposing the dynamical mass constraint improves in- 
deed the quality of the fit, specially for the stellar masses. 
Here we have an other way to look at this, as the evidence 
is a quantitative measure of the aforementioned fit improve- 
ment. The inclusion of the metallicity prior has a similar 
effect of increasing the evidence for the whole sample and 
for the sub-sample of EB stars (both restricted only to stars 
with [Fe/H] measurements). Nevertheless, it does not make 
a big difference for the AS and DK subgroups for which the 
best evidence is still reached in case 2) even when restricting 
only to stars for which [Fe/H] measurements are available. It 
is worth reminding that for most of the DK stars the [Fe/H] 
values we used are the average values for the Taurus- Auriga 
star-forming complex. 

The above analysis suggests that there are still some 
problems with the current generation of standard PMS mod- 
els. However the significance of the disagreement between 
theory and observations is different depending on the sub- 
set of objects considered - for example, we remind that for 
the subset of EBs the average mismatch between observed 
and predicted masses is lower than 10%. Moreover, as it is 
clear from the tables of Appendix|3 for many of the objects 
of the sample the single-star evidence may be the largest for 
other values of the H meta-parameters. Nevertheless this 
global test and the analysis of Sect. |6] both hint to the fact 
that a threefold effort is probably needed to a) improve the 
quality of the data specially assessing the systematic errors, 
b) better constrain the Ay/A.Z, and {Z/X)q values and c) 
improve the physics of stellar models. 



8 APPLICATION TO BINARIES 

In this Section we will analyse the binary systems in our 
data set. For each of them we will check whether the models 
are able to reproduce the coevality of the two stars. In all 
the cases where [Fe/H] measurements are available, we will 
implicitly use the corresponding Gaussian prior on Z in the 
marginalization of the probability distributions. In the four 
remaining cases a fiat prior will be used in the available 
range of metallicities for our models set: Z G [0.007, 0.03]. 

We will mainly make use of the standard H-set previ- 
ously introduced. In the case of severe disagreement between 
the standard models and the observations, we will explore 
the possibility that non-standard H-set might give a better 
agreement with the data. 

Figures [7] and |8] show the HR diagrams with the data 
for each of the 6 EB and 3 AS systems with the best fit- 
ting tracks and isochrones for the standard set of models 
superimposed. The results are also summarized in Table (5] 
Some entries are missing in the table, corresponding to the 
cases in which the confidence intervals are poorly defined. 
This happens when the posterior probability is a very flat 
function and its mode falls outside the confldence interval. 
The best fit masses and ages are obtained by applying two 
different priors in mass, namely a flat and a Gaussian one in 
Figs.[7]and|5] respectively. In the case of EBs, for which stel- 
lar radii are measured, we have used the surface gravity vs. 
effective temperature diagram to compare the models with 
the data. This diagram has the advantage of combining the 
three measured quantities -mass, radius and temperature- 
hence representing the most stringent test for the models. 
For the AS systems we used the HR diagram. We display all 
the best flt models in the HR diagram for homogeneity. 

The best values are obtained after marginalization in 
Z using a Gaussian prior. For displaying purposes only, we 
used isochrones and tracks with a speciflc Z value. The val- 
ues for each system are obtained after transforming the ob- 
served [Fe/H] into Zobs using AY / AZ = 2, Fp = 0.2485 
and {Z/X)q = 0.0181. We then took the closest Z avail- 
able in our models database. In the case where [Fe/H] is 
not available we used Z = 0.0125, the closest to our solar- 
calibrated Z value. The overplotted isochrones correspond 
to the best system composite age, i.e. that obtained by max- 
imizing G'c(-r) = G'p(r) X Gs(r). 

Given the large number and size, the flgures relative to 
the each system's subsection are presented in Appendix iBl 
The upper panels show the marginalized mass distributions 
and the lower ones the age distributions. The left and right 
panels show the results obtained using a flat and a Gaussian 
prior on the mass distribution respectively. 



8.1 RS Cha 

This double-lined E B is located in the Cha cluster 
jMamaiek et al.ll20qol) . The stellar masses and radii are from 
lAlecian et al.l (|2005l ) who refined t he values f rom the pio - 
neering studie s on binary stars by lAndersenI (| 19751 . Il99ll ). 
lAlecian et al.l ([2005) also provide a spectrosco pic measure- 
ment of [Fe/H]. Temperatures are taken from iRibas et all 
(I2OOOI). Surface gravities are simply calculated as \g\ = 
GM/R^. 

Previously thought to be a post-MS system (|jonesl 
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Figure 7. HR diagrams with best fitting tracks and isochrones for the EB and AS binaries of our sample. Masses and ages are obtained 
using a flat mass prior and a Gaussian metallicity prior. 
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Figure 8. HR diagrams with best fitting tracks and isochrones for the EB and AS binaries of our sample. Masses and ages are obtained 
using a Gaussian mass prior and a Gaussian metallicity prior. 
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Table 5. Results from the comparison with the standard set of models. For the missing entries the confidence interval is poorly defined. 



Flat Mass Prior Gaussian Mass Prior 

Mjjyn ^mod Age (stars) Age (system) ^mod Age (stars) Age (system) 

[Mq] [Mq] [Myr] [Myr] [Mq] [Myr] [Myr] 



RS Cha A 


1 


890 


± 


0.010 


1.79 


+0.07 
-0.07 






1.85 


+0.02 
-0.02 


8-oot;;-j^ 


s.ooj 


'0.15 
0.25 


RS Cha B 


1 


870 


± 


0.010 


1.74 


+0.06 
-0.07 


8-70Io;75 




1.82 


+0.01 
-0.03 


_ ^--1-0 30 






RXJ 0529.4 A 


1 


270 


± 


0.010 


1.25 


+0.09 
-0.09 


8.35t?:^^ 


6.25tJ:?» 


1.27 


+0.01 
-0.02 




6.90j 


'1.15 
-0.85 


RXJ 0529.4 B 





930 


± 


0.010 


0.87 


+0.10 

-0.12 


c: 9t- + 1.35 
^•^^-0.70 


0.93 


+0.01 
-0.02 


^■^'-'-0.70 




V1174 Ori A 


1 


009 


± 


0.015 


1.04 


+0.06 
-0.08 


c oc:+0.50 
^•^^-0.40 


7 Qn+0-45 


1.01 


+0.01 
-0.02 


cr oc:+0.50 
O-«O_0.40 


7.40j 


-0.35 
-0.35 


V1174 Ori B 





731 


± 


0.008 


0.42 


+0.08 
-0.07 


„ o^+O 45 




0.73 


+0.00 
-0.02 


o ^ f +0 55 






EK Cep A 


2 


020 


± 


0.010 


1.87 


+0.06 
-0.06 


30.75tf,f 


16.0012-65 


2.02 


+0.00 
-0.02 


26.85l«5f 


18.95 


+ 1.05 
-2.05 


EK Cep B 


1 


124 


± 


0.012 


1.17 


+0.04 
-0.03 


15.80^2.60 


1.13 


+0.01 
-0.01 


18.9011;°^ 




TY CrA A 


3 


160 


± 


0.020 


2.61 


+0.29 
-0.18 




*-^°-0.40 


3.16 


+0.01 
-0.05 




o 7c;+2.65 
'-'•''^-0.20 


TY CrA B 


1 


640 


± 


0.010 


1.52 


+0.24 
-0.35 


3 in+^-55 

■^•^"-0.40 


1.64 


+0.01 
-0.02 


18.901^:05 






ASAS 052821 A 


1 


387 


± 


0.017 


1.54 


+0.08 
-0.09 


J.OU_Q 25 


O.OU_Q 20 


1.39 


+0.01 
-0.02 


o r,c:+0.15 
■^■^'^-0.20 


o 4^+0.10 
■^•™-0.15 


ASAS 052821 B 


1 


331 


± 


0.011 


1.13 


+0.10 
-0.10 


o ^Q+0.15 

O.OU_Q 20 




1.33 


+0.01 
-0.02 


3 60+°-l0 

J.DU_Q 20 






HD 113449 A 





960 


± 


0.087 


0.84 


+0.04 
-0.05 


47.60^^^4^ 




0.86 


+0.04 
-0.04 








HD 113449 B 





557 


± 


0.050 


0.44 


+0.03 
-0.06 






0.48 


+0.02 
-0.04 








NTT 045251 A 


1 


450 


± 


0.190 


1.00 


+0.12 
-0.14 




9 cc+0.65 


1.14 


+0.12 
-0.11 


A 15+3.25 
*-^°-0.90 


3.55j 


-0.85 
-0.50 


NTT 045251 B 





810 


± 


0.090 


0.41 


+0.10 
-0.10 


^•^U_Q 35 


0.65 


+0.07 
-0.08 


3.4o1S;«E; 




HD 98800 Ba 





699 


± 


0.064 


0.51 


+0.25 

-0.02 


n oc+o.oo 


0.85«:?^ 


0.68 


+0.06 
-0.07 


0.851°:?^ 


'^•^'^-o.io 


HD 98800 Bb 





582 


± 


0.051 


0.41 


+0.17 
-0.00 


'J-^3-0.35 




0.56 


+0.05 
-0.06 


I.00l0;;« 







11969 1: 'Ander senlll975l. I1991I ). the X-rays emission reported 
by iMamaick et al. ( 19991 ) clearly points to the PMS nature 
of this system. The two RS Cha components have very sim- 
ilar masses of ~ 1.9M0 and are both close to approaching 
the Zero- Age MS. Recent lite rature estimates for the s ys- 
tem age range fr om Myr (Luhrnan fc Steegh^ 12004 ) to 



9.13 ± 0.12 Myr (|Alecian et al 



200'^ ). 



The outcomes of the comparison between our standard 
set of models and the RS Cha components' gravities and 
temperatures are shown in Fig. mi 

As already noted in Sect. [H] the standard set of mod- 
els slightly underestimates the mass values with /ip = 
1.79lo:o7A^0 and = l.r^to.mMQ for the primary and 
secondary mass respectively. The dynamical masses fall out- 
side these 68% confidence intervals, nevertheless the discrep- 
ancy is quite small -less than 5%- which is a very remark- 
able agreement. The relative precision of the mass estimates 
is strongly increased by the use of the Gaussian prior and 
also the mode of the mass marginal distributions for both 
components are more similar to the observed values with 
MP = 1.83lg:°2A^0 and /is = 1.82t°o:olMQ. 

Even with the slight mass discrepancy, the results on 
the system's age are very robust. The age estimates of the 
two components remarkably agree among each other. From 
the single star's marginal age distribution we obtained the 
combined system age as G(r)RschaA x G(T)RschaB- The 
estimated value for the system age in the case of a flat mass 



prior is rc = 8.50I045 Myr which is narrowed down to 
rc = 8.OOI0 25 Myr when the Gaussian mass prior is im- 
posed. It is worth noticing that the relative precision of the 
combined age is improved with respect to the single stellar 
ages estimates. 



8.2 RXJ 0529.4+0041 A 

The discovery of this double-lined EB located in th e Orio n 
star-forming region was reported by ICovino et al. I (|2000t ). 
The same group refined the system parameters using new 
photometric observations in lCovino et alj (|2004 ). We adopt 
the data from the la tter paper and the Orion [Fe/H] from 
iD'Orazi et al] (|2009l ). 

Comparing the observed gravity and temperatures with 
our standard set of models yields stellar masses in agreement 
with the dynamical measurements with /ip = 1.25lo ogMQ 
and fis = O.STIq 5^2-^0 when a flat mass prior is used. If 
a Gaussian prior is applied, then the precision improves 
by a factor of 10 for the estimated masses with fxp = 



1.27lo.o2^0 and /is = O.93lo.o2^0- On the other hand the 
derived single stellar ages are in slight disagreement. Using a 
Gaussian mass prior, we obtain respectively rp = 8.70l|^'25 
Myr and rs = 5.20lJ'7o Myr. Nevertheless, the two ages 
are both consistent with the composite age of the system, 
i.e. rc = 6.90log5 Myr. The results for the standard set of 
models are summarized in Fig. IB2I 



3+0.01 
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By looking at the Bayes factors of Table IA4I it appears 
that most of the age discrepancy might be ascribed to the 
poorer fit of the secondary star. For the primary star the 
standard set of models provides the second highest Bayes 
Factor value, 18.34, the highest being just 18.35 for the set 
with AF/AZ = 2,q: = 1.68 but Fp = 0.23. For the sec- 
ondary the standard set provides a Bayes Factor that is 
~ 1.5 times smaller than the one giving the largest evi- 
dence, i.e. Ay/AZ = 2, a = 1.20 but Fp = 0.2485. This is 
not enough to state that the latter set gives a significantly 
better agreement with the data, but using the latter set of 
meta-parameters yields an age of 6.30lj gg Myr, which is 
in agreement with the primary star's age within the uncer- 
tainty interval. 

The fact that the secondary star is better fitted 
by cooler model s (i.e. models with l ower a) was al- 
ready reported bv l D'Antona et al] l|2000l ') and confirmed in 
ICovino et all \2Qm . 



8.3 V1174 Ori 



This double-lined EB was discovered by IStassun et al.l 
l|2004ll . We adopt stellar parameters fr om this paper and th e 
average [Fe/H] abundances for Orion jP'Orazi et al.|[2009l ). 

As in the case of RXJ 0529.4-1-0041 A, the primary 
star of V1174 Ori is moving away from the Hayashi track, 
while the secondary is still fully convective. Also in this 
case the models show some difficulty in reproducing the 
secondary observables. The standard set of models -with 
a flat mass prior- predicts a secondary mass of ~ 0.42Mq, 
much smaller than the dynamical mass (~ 0.7Mq). The pri- 
mary mass is instead well recovered with /ip = I.mIq o8-^0 • 
The situation for the secondary does not improve much even 
when using the coldest set of models available. 

It has been noted (see e.g. iHillenbrand fc White! 120041 ') 
that one of the problems in estimating the effective temper- 
atures for PMS stars from the observed spectral type is the 
adoption of temper ature scales that are calibrated on MS 
stars. For example lStassun et al.l (2004) use the temperatur e 
scale for dwarf stars by Schmidt-Kaler in lAUer et al.l l|l982l ). 
The same authors show how stellar models are not able to 
reproduce luminosities and temperatures for the secondary 
star and attribute the discrepancy to the non-adequacy of 
the dwarf spectral type to T^s conversion when applied to 
PMS stars. 

At a given spectral type PMS stars are in general hotter 
than the corresponding MS stars. Their surface gravities in- 
dicate that a temperature scale intermediate between dwarfs 
and giants should be adopted. IHillenbrand fc White! (|2004! ) 
suggest that temperature corrections as high as 100 K could 
be necessary to compensate for the temperature underesti- 
mates. We explored this possibility in the case of V1174 Ori, 
by artificially increasing the temperature of the primary by 
100 K and keeping the effective temperature ratio between 
primary and secondary constant. The reason why we can not 
simply increase the effective temperature of the secondary, 
which is the main responsible of the d isagreement with the 
models, is that in !Stassun et al.l (|2004! ') this quantity is not 
directly and independently measured, but it follows from the 
determination of the primary effective temperature from the 
spectral type and the temperature ratio from the light curve. 



Hence, if any offset is present, it should be reflected in both 
components. 

The results of the comparison of the modified observ- 
ables with the standard set of models are shown in Fig. IB4I 
The situation is only slightly improved compared to Fig. lB3l 
The gain in the secondary mass estimate is that now the 
best mass is ~ 0.49Mq, not yet enough to be in agreement 
with the dynamical mass. However the primary mass is still 
recovered within the uncertainty interval. It is clear that a 
change in Tcff has a larger impact on the inferred stellar 
mass when a star is still in the vertical Hayashi track than 
when it is located on the almost horizontal Heyney track. 

Even with the small improvement achievable by increas- 
ing the estimated Tcff, V1174 Ori remains a challenge for 
stellar evolution theory. However parallel observational ef- 
forts are required to assess the issues related to the effective 
temperatures determinations. 

8.4 EK Cep 

This system is kn own to be an EB for more than 50 years 
(jStrohmeiei! Il959! ) . Quite some time after its discovery it 
was recognized to host a ~ I.IMq PMS star, together 
with a ~ 2.0 M(7^ primary already on its main sequence 
ijPoppei! 1 1987! ). We adopted stellar parameters from this 
paper and the spectro scopic determination of [Fe/H] by 
llMartin fc Rebold (|l993l ). 

With our standard set of models the stellar masses for 
EK Cep are not recovered within the 68% confidence in- 
terval (see Fig. IB5p . If we use a fiat mass prior we ob- 
tain /ip = l.d,7tollMQ and /xp = 1.17i;;:^3M0, while 
the measured dynamical values are 2.020 ± O.OIOMq and 
1.124 ± 0.012Mq for primary and secondary, respectively. 
Hence the primary mass is slightly underestimated and the 
secondary mass slightly overestimated. It has to be noted 
that the absolute difference between model-predicted and 
dynamical masses are of the order of 7% and 4% hence quite 
small. Nevertheless, the discrepancy are significant accord- 
ing to our deflnition of the confidence interval. 

Given the Bayes Factors of Table 1X41 we have compared 
the data with models from the class with a — 1.20, Yp — 
0.23, and AY/AZ = 2, i.e. the one that gives the largest 
composite evidence for the system. The agreement is still 
not satisfactory. As shown in Fig. IB6I a further improvement 
of the fit is achieved by assuming a higher metallicity (i.e. 
Z = 0.01 93 rather than Z = 0157) as ff the old {Z/X)q = 
0.0231 bv lGrev esse fc S auva il (ll998!) was used instead of the 



recent one by Asplun d et al.l 1 2009! ). The inferred masses 



+0.14 



M(7 



are in this case np = 1.97lQ QgMQ and np = 1.12_q q^im^ 
This r esults are in very good agreement with those of lClaretl 
||2006! ). who using similar values, namely Z = 0.0175 and 
a = 1.3, is able to reproduce the system observables. 

This test shows how the success of a set of models in 
reproducing the observations might be severely affected by 
the current uncertainties on the meta-parameters. Paradox- 
ically, in the CEise of EK Cep system models calculated with 
out-of-date meta-parameters seem to give a better agree- 
ment with the data than the state-of-the-art ones. 

One interesting thing that this system shows about our 
method is the power of the combined system age marginal 
distribution, G'(r)EKCopA x G'(t)ekCopB. Since the primary 
star is already on the Main Sequence, its evolution is very 



Testing pre-main sequence models: the power of a Bayesian approach. 



19 



slow, resulting in a very flat G(t) and, consequently a very 
poor precision in the age determination. Nevertheless the 
age of the system is very well determined -given the choice 
of the model class. As a consequence also the primary star 
has a very precise age determination, which is very valuable 
for MS stars. In the case of the standard set of model and 



Gaussian mass prior the system age is r = 18.953 



Myr, 



while for the non-standard set used in this Section, again 
with a Gaussian mass prior, we obtain t 



26.55l?:«^Myr. 



8.5 TY CrA 



This double lined EB is part of a hierarchic al system with 
three or possibly four stellar components (see lChauvin et al] 
I2OO3I). The fundamenta l parameters we adopted are taken 
imfc 



from iCasev et all (jl998') . For this particular system we could 
not find any spectroscopic determination of [Fe/H]. There- 
fore, instead of applying a Gaussian prior on the metallicity, 
the marginalization over Z was made using a flat prior with 
Z G [0.007,0.03], i.e. the range of metallicities available in 
our models grid. 

As in the case of EK Cep, the primary star is already on 
the Main Sequence, while the slower evolving secondary is 
still on its Hayashi track. Similarly to the EK Cep case, our 
standard set of models is able to reproduce the secondary 
mass quite well, while the primary mass is once again un- 
derestimated. The values we obtain when a flat mass prior 
is used are = 2.611q'i8-^0 ''■'^d Ms ~ 1-52 



M(7i while 



the dynamical masses are estimated to be 3.16 ± 0.02AfQ 
and 1.64 ± O.OIMq for the primary and secondary, respec- 
tively. The low relative precision of these model predictions 
-compared e.g to the case of EK Cep- are mainly due to the 
larger uncertainties on the effective temperatures and radii 
for the TY CrA system (see Table H]). 

Also in this case, we tried to see whether the agree- 
ment between data and models might be improved by using 
the coldest set of models with a = 1.20, Ip = 0.23 and 
AY/ AZ = 2. As it is possible to see by comparing Figs. 
[B7]and[B8l there is a slight improvement in the primary 
mass determination, without losing the good agreement for 
the secondary mass. The results obtained with a flat mass 
prior are fip = 2.691EJ:?^Mq and fis = lAdtHtMQ. The 
primary mass is still underestimated by about 15%, not a 
too bad result -in absolute terms- but further investigation 
is needed to explain this partial disagreement. 

TY CrA A is on the MS, therefore its age is not well 
constrained. Nevertheless the mode of the age distribution 
for the primary is very close to the mode of the secondary, 
which has a better constrained age determination. From the 
composite age distribution G(r)TYCrAA x G(r)TYCrAB in 
the case of the standard set and applying a Gaussian mass 
prior we obtain r = 3.75to |^ Myr. For the coldest set of 
models and still applying a Gaussian mass, prior we obtain 
a slightly older age of r = 5.20t.Q'jo My r. These age values 
are slightly older than the age found bv lCasev" et al.l (|l998l ') 
who roughly estimate a system age of ~ 3 Myr. Also these 
authors show that the models have problems in consistently 
predicting the stellar observable for both components. While 
being able to reproduce the observed secondary properties, 
they also find that models overestimate the effective temper- 
ature of the primary star. This is equivalent to our finding of 
an underestimated stellar mass for the given Tcfi and \ogg. 



8.6 ASAS J052821+0338.5 

This is the most recently discovered double- lined EB in our 
data set. The stellar parameters are from IStempels et al.l 
(2008i). The two stars have very similar masses, and the 
slightly more massive primary (1.387 ± 0.017Mq) is just 
moving away from its Hayashi track while the secondary 
(1.331 ±O.O11M0) is located just before the end of the fully 
convective phase. The [Fe/H] value we used is —0.15 ± 0.2 
dex, i.e. the average of the quoted values for the primary 
(-0.2 ± 0.2 dex) and the secondary (-0.1 ± 0.2 dex) in 
IStempels et al.l (|2008l ). 

The predictions of our standard set of models slightly 
differ from the measured masses. The values we obtain 
when a flat mass prior is used are /xp = 1.54j;Q jjMQ and 
= l-13lg jQMQ. Therefore the primary mass is overes- 
timate d by 11% and the se condary is underestimated by 
~ 15%. IStempe¥ et al.l l|2008f ) provide a double solution for 
the system parameters depending whether stellar spots are 
taken into account in the light-curve analysis (as in the case 
reported in Table |4| or not. We applied our method using 
also the measurements for the latter case and the results are 
shown in Fig. lBlOl Both of the predicted masses are in better 
agreement with the data in this case with /ip — 1.53to io-^0 
and fis = 1.241q JgMQ. Hence for the primary mass the sit- 
uation is slightly better with an overestimate of 10%, while 
the situation is much improved for the secondary which is 
now predicted to be 7% less massive than the observed value. 

It is clear that the detailed modeling of the light curve 
plays an important role in determining stellar properties 
and, as a consequence, in constraining the models' predic- 
tions. Once again the models give an overall satisfactory 
agreement, being 10 or even 15% still a quite good error in 
stellar masses predictions. Nevertheless more work is needed 
to explain these differences. 

The results for the age of the system are more robust 
with the two components having ages in good agreement. In 
the case where the light curve solution including star spots 
is used, we obtain a system age -using a Gaussian mass 
prior- of 3.45^0 Myr. When using the light-curve solu- 
tion without star spots, we instead obtain 3.651q'2o M!yr. In 
both cases the age is much younger than what is found by 
IStempels et al. (I2OO8I) who, using solar metallicity models 
bv lBaraffe et all l| 19981 ). found an age of ^ 10 Myr for the 
system. We used the [Fe/H]~ —0.15 dex quoted in the same 
paper to derive our Gaussian Z prior. This value is slightly 
sub-solar, hence p art of the difference be tween our age esti- 
mate and that by IStempels et al] (|2008l ) could be ascribed 
to that. Still it is quite hard, even using our solar metallicity 
models, to reproduce a ~ 1.4Mq star close to the base of 
its Hayashi track at such an old age like 10 Myr. 



8.7 HD 113449 

This system is an AS binary whose orbital param eters have 
been recently estimated by ICusano et all i|20ld ). Here we 
use slightly different parameters (yet unpublished) kindly 
provided by the same grou p after more accurate an alysis of 
the data and the [Fe/H] b ylPaulson fc Yeldal (120061 '). 
As already noted by ICusano et al. 

I l|2010i r there is a 

slight disagreement between dynamical and inferred masses 
by several sets of stellar models. Also the masses pre- 
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dieted by our standard set of stellar tracks are slightly 
underestimated. The primary mass is found to be /xp = 
QMtoltMQ while the secondary is /^s = O.44j:5;;^^M0 
(see also Fig. IBTT]| . From Table [5] we can see the primary 
dynamical mass of 0.960 ± O.OSTMq is still consistently 
recovered while for the secondary the dynamical mass of 
0.557 ± 0.050AfQ is outside the 68% confidence interval. 

To see whether the discrepancy could be reduced, we 
used the coldest set of models and, in addition, we de- 
rived the Z value usin g the spectroscopic JFe/H ] and the 
{ZIX)q = 0.0231 bv ICrevesse fc Sauvall l|l998h . The re- 
sults, displayed in Fig. IB12I show a better agreement 
with the observations. The improvement is not substantial, 
though, and the predicted stellar masses are in this case 
/^p = O.89lH5^0 and = QA^t°of^MQ. 

The stellar ages in this particular case are not very well 
determined. The two stars are indeed very close to their MS 
position, which makes age determination very hard. Never- 
theless the primary shows a small peak in its G{t) distri- 
bution at an age of r ~ 50 Myr (for both the standard 
and non-standard set of models). The secondary instead 
does not show any peak in the stellar age, with a very flat 
G(r) slightly increasing towards the edge of our models age- 
interval (100 Myr). The system age is poorly defined as well. 

8.8 NTT 045251+3016 

The discovery of this AS binary was first reported by 
ISteflten et all l|200ll ') from which we adopted the stellar pa- 
rameters. In this case no spectroscopic [Fe/H] is available. 
This system is quite young and both the primary and the 
secondary are found in their fully convective phase a long 
the Hayashi track. As pointed bv ISteffen et al.l ()200ll ). all 
the stellar models adopted by them predict too low masses 
for both components. The set of models that gives the best 
agreem ent with observations is the one by iBaraffe et al.l 
l|l998h when a mixing length parameter a = 1.0 is adopted. 
This is not surprising given that models with a lower a, be- 
ing intrinsically colder, predict larger masses for given ob- 
served luminosities and temperatures. 

With our standard set of models and using a flat mass 
prior, the two masses are found to be /ip = l.OOlg j^Mq 
and /is = 0.411q ;[qMq, severely lower than the dynamical 
masses by ~ 30% and ~ 50 % respectively (see Fig. IB13)) . 
From Tables lAll and IA2I it is possible to see that for both 
of the components of NTT 045251+3016 the set of models 
that provides the largest Bayes Factor is, once again, the 
coldest set available (Yp = 0.23, AY/AZ = 2, q = 1.2). 
When using this particular set and a flat mass prior, we 
obtain slightly larger masses of /ip = 1.13lQ ;t|MQ and 
/is = 0.50lo i2AfQ (see Fig. IB14|I . The improvement is 
not enough to obtain an agreement between predicted and 
observed mass for the secondary, while the primary mass, 
though still underestimated, is in agreement within the er- 
rors. 

We report that using the BASE software (courtesy of 
Tim Schulze-Hartung, in prep) for analysing the system's as- 
trometric measurements and radial velocities, slightly lower 
masses are predicted. The primary mass is found to be 
/ip = 1.383 ± 0.220Mq (-4.60%) while the secondary is 
/is = 0.766 + 0.089 (-5.41%). 

Even with this latter improvements, there is still a larger 



disagreement in the predicted vs. dynamical mass for NTT 
045251+3016 than what we found for the EBs cases or even 
the other two AS binaries. This suggests that part of the 
problem might reside in observations as well, and we already 
noted in Sect. [5] the peculiar location of this star in the HR 
diagram given its measured mass. Apart from a theoretical 
effort, which is certainly needed, this system demands atten- 
tion also from the observational side to exclude, e.g., higher 
order multiplicity that to-date interferometric observations 
are not capable to resolve. 

Concerning the system composite age, the values we 
obtain when using a Gaussian mass prior are r = 3.55^0 50 
Myr and 4.65^0 65 Myr for the case of standard and coldest 
set of models, respectively. 

8.9 HD 98800 B 

This AS binary is part of a quadruple system. iBoden et al.l 
('2005') reported preliminary visual and physical orbit for the 
binary subsystem. They derived the components' masses of 
0.699 + O.O64M0 and 0.582 + O.O51M0 for the primary and 
the s econdary, respectiv ely. We adopted the [Fe/H] value 
from iLaskar et all l|2009l '). Both the components of the sys- 
tem are very young and located at the beginning of their 
Hayashi track. 

For this system the standard set of models provides 
a good fit to both components (see Fig. IB15|I . The pre- 
dicted mass values in the case of a flat mass prior are 
/ip = O.511!5;o|M0 and /is = O.4lto.Jo-^0- The location 
of the two stars in the HR diagram makes their marginal 
distribution extremely asymmetric. As a consequence, the 
best values are located quite close to (or exactly at) the 
boundary of the confidence intervals. This is the reason why 
the quoted lower error for the secondary mass is zero. The 
best values for both the primary and secondary are slightly 
smaller than the dynamical mass values but in this case 
there is consistency within the 68% confidence intervals. 

The inferred ages for the two components are very sim- 
ilar and when a Gaussian mass prior is adopted we obtain 
rp = 0.85to;?o Myr and rs = l.OOlHs Myr. Also in this 
case the marginal distributions are quite asymmetric. The 
composite system age is found to be rc = 0.901q'5o Myr, 
very young indeed. 



9 THE STARS IN TAURUS-AURIGA 

In this section we will present the results for the 9 stars 
found in the Taurus- Auriga star-forming region and whose 
masses are derived using disk kinematics (see Sect.[5)|. For 
all of these stars, the Gaussian metallicity prior is applied 
using the average value of [Fe/H] = —0.01 + 0.05 dex for the 
region (D'Orazi et al. 2011). 

This sample exactly coincides with the DK+DKS sam- 
ple of Sect. [5] As it is possible to see in the rightmost panels 
of Fig. |6l most of the DK+DKS stars have strongly under- 
estimated values of the mass. We discussed some possible 
reasons for this discrepancy in Sect.|6] 

We used the standard set of model to derive the ages of 
the Taurus- Auriga stars. The predicted values when a Gaus- 
sian mass prior is used are reported in Table [6l In addition 
to single stellar ages we also computed the composite age 



Testing pre-main sequence models: the power of a Bayesian approach. 



21 



Table 6. Derived ages for Taurus-Auriga DK stars. 



Name 


Age[Myr] 


Rel. 


Prec. 


UZ Tau Ea 


2-65Io.45 


0, 


,373 


UZ Tau Eb 




0, 


,592 


UZ Tau E" 


l.SO_Q 25 


0, 


,199 


DL Tau 


1 20+"-^° 


0, 


,195 


DM Tau 


■^•^■^-0.50 


0, 


,276 


CY Tau 


1 85+"-='° 


0, 


,338 


BP Tau 


1 7 oe;+S.85 
^'■^3-3.35 


0, 


,370 


GM Aur 




0, 


,095 


MWC 480 


10.30*' 


0, 


,957 


LkCa 15 


^•^^-0.75 


0, 


,371 



" Value for the composite system age 
^ Uncertainty interval poorly defined 

for the DKS system UZ Tau E. In general these age deter- 
minations have a worst precision when compared to the EB 
and AS sample. Moreover for the star MWC 480 the peak of 
the distribution is outside the 68% confidence interval. This 
is because the G(r) is very fiat and the corresponding peak 
is just barely visible. This peak is located in the area cor- 
responding to the leftmost 16% probability that is excluded 
according to our definition. This is not strange given our 
definition of the confidence interval, it is just an indication 
that the age of this system is very poorly defined. 

As we mentioned in Sect.[S]the stars UZ Tau Ea, BP Tau 
and MCW 480 have a peculiar location in the HR diagram. 
We excluded the latter two stars, while we kept UZ Tau Ea 
in the sample to obtain an average age of 2.1 ± 1.3 Myr 
where the quoted uncertainty is the standard deviation of 
the ages of the remaining stars. For both the stars in the UZ 
Tau E system we considered the composite age as the best 
age estimator. This average age is in very good agreement 
with the estimated age for the Taurus- Auriga star for ming 
region of 1-2 Myr (see e.g iKraus fc Hillenbrancill2009l . and 
references therein). 



10 SUMMARY AND CONCLUSIONS 

The importance of a stringent test of PMS models against 
stars with accurately known parameters (i.e. mass, lumi- 
nosity, radius, effective temperature, [Fe/H]) can hardly be 
overestimated, as these models represent the main tool to 
derive masses and ages of stars observed in star forming re- 
gions and young stellar clusters. Consequently the inferred 
star formation histories and mass functions of young stellar 
groups strongly depend on the answers provided by stellar 
evolutionary codes. 

In order to constrain PMS models we relied on a data 
set containing 25 PMS stars of measured mass (plus 2 MS 
companions in binary systems). This is the full up-to-date 
sample of known PMS stars with dynamical mass measure- 
ments in the range 0.2-3.0 Mq. Among them 10 PMS ob- 



jects belong to double-lined eclipsing binary systems and 6 
to astrometric and spectroscopic binaries; the remaining 9 
objects are stars whose masses are derived using the mea- 
sured orbital velocity of their circumstellar disks. 

The main novelties of the paper are both the approach 
followed for comparing theory with observations and the set 
of PMS models used in the comparison. Regarding the for- 
mer, we applied for the first time to the whole sample of 
PMS stars a very general Bayesian method. This approach 
allows a full exploitation of the available information about 
the observed objects which is included in the form of prior 
probability distributions. In addition it provides robust un- 
certainties for the inferred quantities. 

The mode ls are extracted fro m the very recent Pisa 
PMS database l|TogneUi et al.ll201lh . They include the state- 
of-art input physics and are available for a large and very 
fine grid of metallicities, masses and ages and for different 
primordial helium abundances, Yp, helium-to-metals enrich- 
ment ratios, AY/ AZ, and mixing-length parameter, a, val- 
ues. 

We checked the robustness and accuracy of the method 
in recovering stellar ages and masses against simulated bi- 
nary data sets. One interesting result is that even synthetic 
binary stars -coeval by construction- might mimic non co- 
evality as a consequence of the random uncertainty in the ef- 
fective temperature, radius, and luminosity. The actual frac- 
tion of fake non coeval recoveries strongly depends on the 
sample characteristics, since we demonstrated that the abil- 
ity of recovering the simulated masses and ages is a complex 
function of the actual position of the star in the HR dia- 
gram. This fraction can be as large as 95% for systems with 
one component close to its main sequence position. 

This suggests that the inability to fit both components 
of a binary system with a single isochrone does not neces- 
sarily imply that the two stars are not coeval or that the 
models present some deficiency. We also showed that, even 
in the simulated systems that are recovered as non coeval, 
the inferred system age obtained using the composite age 
distribution is in very good agreement with the simulated 
age. 

When the real data are used, the Pisa PMS models 
show an overall agreement with the observations. With the 
exception of V1174 Ori B, the masses of EB stars are well 
recovered within 10%. The agreement progressively worsens 
for AS binaries and DK stars, but also the observational 
uncertainties become more severe for the latter objects. 

A slightly worse situation is observed for stellar ages 
compared to stellar masses. Within our sample 6 binary 
systems are present for which ages can be derived for both 
components. The systems EK Cep and TY CrA both have a 
primary on the MS, hence their primaries' ages are not very 
accurate and we exclude them from the following discussion. 
HD 113449 is excluded as well, since the age of the secondary 
is also not well determined. Of the remaining 6 systems, 4 
are predicted to be coeval -within the uncertainties- by our 
models. The systems RXJ 0529.4-f0041 A and V1174 Ori 
are instead not consistent with coevality. 

We have shown by simulations that the probability of 
deriving non coeval ages for two stars in a binary system 
-because of observational uncertainties- is a complex func- 
tion of the position of stars in the HR diagram. Nevertheless 
we have also shown that, for the subset of simulated system 
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for which both components are stiU far from their main se- 
quence position, this probabiUty is very low, always below 
10% and as small as 0% for most simulated cases. The ob- 
served non coevality fraction of 30% is higher than what 
expected, given the typical quoted observational errors and 
the location in the HR diagram of the stars of the 6 afore- 
mentioned systems. With the present low number statistics, 
it is hard to draw robust conclusions, nevertheless these re- 
sults might suggest that PMS models still have some prob- 
lems in correctly predicting stellar ages. On the other hand 
the hypothesis that some systems might really be non co- 
eval can not be ruled out given the current uncertainties on 
the models. An other possibility is that observational errors, 
specially the systematic ones are underestimated, leading to 
a spurious non coevality of the stars in RXJ 0529.44-0041 A 
and V1174 Ori. 

With our Bayesian approach it is possible to evaluate 
the probability for different sets of models, i.e. the mod- 
els' evidence. We analysed the entire data set using several 
classes of models computed with different Ip, AY/AZ, and 
a values. We calculated the evidence for each star using 9 
different meta-parameters configurations. Furthermore, four 
combinations of the prior distributions for mass and metal- 
licity have been used for each meta-parameter choice, for a 
total of 36 classes. We found that adopting a Gaussian rather 
than a flat mass prior significantly improves the composite 
evidence for the full data set; the same effect, but to a lesser 
degree, is obtained imposing a Gaussian metallicity prior, 
mainly for EBs. 

Although our standard set of models shows a reasonable 
general agreement with the data, predicting mass values al- 
most always within 20% of the dynamical ones -and in sev- 
eral case even within 5%-, the general trend suggests that 
standard models tend to under estimate the stellar mass, 
confirming previous results (see iMathieu et al.l |2007| . and 
references therein). As a consequence, the largest composite 
evidence is obtained with our coldest set of models, i.e. with 
the mixing-length parameter q=1.2 and the lowest helium 
abundance at fixed metallicity. 

Given that the discrepancy between theory and obser- 
vations increases going from the most precise data set of EBs 
to the others, we point out that a twofold effort is needed 
to achieve a better agreement. From the theoretical point of 
view a better understanding on the treatment of superadi- 
abatic convection and a better characterization of the mod- 
els' meta-parameters is desirable. Moreover, one of the open 
questions of PMS modeling concerns the role played by ac- 
cretion in affecting the observable properties of very young 
stars, i.e. their effective temperatures and luminosities. Dif- 
ferent authors disagree on the impact accretion can have 
on the s tellar ages which are inferred using non- accreting 
models (|Baraffe et al.ll2009l : ITlosokawa et aLlbOllT ). Unfor- 
tunately, a fully consistent treatment of accretion for stars 
spanning a wide range of parameters is still lacking, there- 
fore preventing a proper and quantitative evaluation of the 
impact of accretion on the inferred stellar properties. 

From the observational side the significance of such a 
comparison could be improved in the future by a larger sam- 
ple of well studied and characterized PMS stars and by a 
better control on the systematic errors affecting AS and DK 
stars' measurements. 
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Table Al. Ba\-('s fac tor \-ahi('s foi' the case of flat prior on the iriass and flat prior on the nietallicity 



a = 1.20 a = 1.68 a = 1.90 



a = 1.20 a = 1.68 a = 1.90 



AY/AZ -- 
AY/AZ -- 
AY/AZ -- 



: 2, Yp : 
: 2, Yp : 
: 5, Yp : 



0.23 

: 0.2485 
: 0.2485 



RS Cha A 



1.252 
1.004 
0.556 



1.252 
1.000 
0.561 



1.253 
1.009 
0.559 



RS Cha B 



1.054 
0.975 
0.787 



1.077 
1.000 
0.802 



1.111 
1.033 
0.828 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
2, Yp : 
: 5, Yp : 



RXJ Aa 



RXJ Ab 



0.23 
: 0.2485 
: 0.2485 



1.002 
0.955 
0.833 



1.042 
1.000 
0.889 



1.035 
0.996 
0.889 



1.377 
1.292 
1.038 



1.101 
1.000 
0.716 



1.003 
0.900 
0.620 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



V1174 Ori A 



V1174 Ori B 



0.23 
: 0.2485 
: 0.2485 



1.108 
1.073 
0.960 



1.051 
1.000 
0.824 



1.009 
0.948 
0.741 



1.589 
1.365 
0.821 



1.183 
1.000 
0.579 



1.086 
0.914 
0.525 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp 
2, Yp 

5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



EK Cep A 



0.950 
0.995 
1.451 



0.955 
1.000 
1.457 



0.956 
1.005 
1.460 



EK Cep B 



5.684 
4.447 
1.386 



1.181 
1.000 
0.712 



0.946 
0.866 
0.712 



AY/AZ : 
AY/AZ : 
AY/AZ 



2, Yp : 

2. Yp 



: 0.23 

0.2185 
0.2f<sr, 



TY CrA A 



0.962 

0.999 
l.lfr, 



0.963 
1.000 
1.1 f7 



0.964 
1.001 

1.117 



TY CrA B 



1.166 
1.111 
0.9.^)7 



1.066 
1.000 
0.81.-) 



0.977 
0.909 
0.720 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



ASAS A 



ASAS B 



: 0.23 
: 0.2485 
: 0.2485 



0.922 
0.884 
0.784 



1.036 
1.000 
0.906 



1.052 
1.017 
0.913 



1.251 

1.208 
1.071 



1.083 
1.000 
0.745 



0.943 

0.848 
0.580 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



HD 113449 A 



0.23 
: 0.2485 
: 0.2485 



0.569 
0.667 
1.003 



0.971 
1.000 
0.837 



HD 113449 B 



1.019 
0.998 
0.670 



0.665 
0.806 
1.272 



0.832 
1.000 
1.589 



0.883 
1.060 
1.686 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
: 2, Yp : 
: 5, Yp : 



0.23 
: 0.2485 
: 0.2485 



NTT A 



1.422 
1.340 
1.101 



1.098 
1.000 
0.723 



0.976 
0.877 
0.607 



NTT B 



1.776 
1.529 
0.946 



1.173 
1.000 
0.612 



1.043 
0.887 
0.541 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
2, Yp : 
: 5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



HD 98800 Ba 



3.162 
2.816 
1.937 



1.120 
1.000 
0.754 



0.841 
0.756 
0.600 



HD 98800 Bb 



2.719 
2.322 
1.563 



1.148 
1.000 
0.911 



0.903 
0.811 
0.632 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



UZ Tau Ea 



UZ Tau Eb 



0.23 
: 0.2485 
: 0.2485 



1.731 
1.540 
1.202 



1.126 
1.000 
0.795 



0.968 
0.846 
0.623 



2.795 
2.228 
0.897 



1.363 
1.000 
0.229 



0.862 
0.614 
0.123 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp 
2, Yp 

5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



DL Tau 



2.465 
2.147 
1.400 



1.139 
1.000 
0.735 



0.924 
0.814 
0.616 



DM Tau 



1.784 
1.529 
0.929 



1.182 
1.000 
0.589 



1.051 
0.886 
0.518 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



CY Tau 



1.551 
1.371 
0.996 



1.131 
1.000 
0.735 



1.025 
0.902 
0.646 



BP Tau 



1.103 
1.042 
0.890 



1.058 
1.000 
0.848 



1.045 
0.987 
0.832 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 

2, Yp : 
5, Yp : 



GM Aur 



: 0.23 

: 0.2485 
: 0.2485 



7.963 

6.900 
5.211 



MWC 480 



1.329 

1.000 
0.244 



0.287 
0.173 
0.023 



1.118 

0.991 
0.564 



1.128 

1.000 
0.570 



1.131 

1.002 
0.571 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 

2, Yp : 
5, Yp : 



0.23 
: 0.2485 
: 0.2485 



LkCa 15 



1.259 
1.202 
1.029 



1.080 
1.000 
0.754 



0.992 
0.904 
0.646 



Testing pre-main sequence models: the power of a Bayesian approach. 

Table A2. BaycK fac tor \-ahi('s foi' the case with Gaiissian prior on the mass and flat prior on the niotallicity 







a = 1.20 


a = 1.68 


a = 1.90 




a = 1.20 


a = 1.68 


a = 1.90 


AY/AZ = 2, Yp = 0.23 


Ho Una A 


6.968 


7.171 


7.039 


tib Ona D 


17.54 


17.92 


18.36 


AY/AZ = 2, Yp = 0.2485 




7.123 


7.150 


7.212 




12.47 


12.64 


12.88 


AY/AZ = 5, Yp = 0.2485 




0.147 


0.142 


0.143 




0.003 


0.003 


0.003 


AY/AZ = 2, Yp = 0.23 


RXJ A 


11.72 


16.17 


13.57 


RXJ B 


11.71 


13.11 


12.30 


AY/AZ = 2, Yp = 0.2485 




8.582 


16.72 


14.98 




12.81 


11.93 


10.14 


AY/AZ = 5, Yp = 0.2485 




2.719 


14.44 


15.56 




12.08 


3.470 


1.650 


AY/AZ = 2, Yp = 0.23 


V1174 Un A 


2.844 


7.566 


9.745 


V1174 (Jri B 


6.218 


1.236 


0.658 


AY/AZ = 2, Ip = 0.2485 




4.949 


10.52 


12.47 




3.032 


0.358 


0.158 


AY/AZ = 5,Yp = 0.2485 




16.15 


12.65 


6.966 




0.024 


< 10-3 


< 10-3 


AY/AZ = 2, Yp = 0.23 


EK Cep A 


9.530 


9.597 


9.580 


EK Cep B 


18.02 


17.25 


9.205 


AY/AZ = 2, Yp = 0.2485 




2.259 


2.254 


2.256 




17.48 


23.61 


14.84 


AY/AZ = 5, Yp = 0.2485 




0.001 


0.001 


0.001 




28.67 


14.20 


21.99 


AY/AZ = 2, Yp = 0.23 


1 Y Or A A 


0.409 


0.409 


0.409 


1 Y UrA D 


13.52 


6.466 


4.870 


AY/AZ = 2. Yp = 0.2485 




0.145 


0.145 


0.115 




13.13 


6.728 


4.847 


AY/AZ = .'). 5p = 0.2 l.s.') 




O.Ofl.') 


0.0 or) 


0.005 




')!')(} 


().8():-) 


1.718 


AY/AZ = 2, Yp = 0.23 


A O A O A 

AbAb A 


3.446 


0.909 


2.070 


AoAo B 


8.715 


12.54 


7.418 


AY/AZ = 2, Yp = 0.2485 




9.576 


2.131 


3.682 




11.10 


9.666 


3.843 


AY/AZ = 5, Yp = 0.2485 




29.33 


7.729 


12.42 




21.51 


0.383 


0.023 


AY/AZ = 2, Yp = 0.23 


tiu 113449 A 


5.754 


10.44 


10.94 


T TTA ~i ~i ^ A A r\ T~> 

HD 113449 B 


2.548 


3.270 


3.497 


AY/AZ = 2, Yp = 0.2485 




5.856 


9.333 


9.159 




2.548 


3.252 


3.476 


AY/AZ = 5, Yp = 0.2485 




4.563 


2.854 


2.015 




2.559 


3.129 


3.291 


AY/AZ = 2, Yp = 0.23 


JN i i A 


3.589 


1.323 


0.830 


■\Trprp T~) 
JN 1 i B 


2.417 


0.372 


0.192 


AY/AZ = 2, Yp = 0.2485 




2.837 


0.880 


0.516 




1.391 


0.166 


0.078 


AY/AZ = 5, Yp = 0.2485 




1.104 


0.173 


0.079 




0.116 


0.005 


0.002 


AY/AZ = 2, Yp = 0.23 


HD 98800 Ba 


2.327 


6.120 


5.958 


HD 98800 Bb 


4.910 


7.595 


7.702 


AY/AZ = 2, Yp = 0.2485 




2.909 


6.173 


5.419 




4.354 


7.662 


6.594 


AY/AZ = 5, Yp = 0.2485 




4.946 


5.024 


3.010 




7.062 


6.963 


3.018 


AY/AZ = 2, Yp = 0.23 


UZ iau hid. 


0.077 


< 10-3 


< 10-3 


UZ lau ejD 


46.45 


16.35 


9.575 


AY/AZ = 2, Yp = 0.2485 




0.030 


< 10~3 


< 10-3 




33.43 


9.941 


5.498 


AY/AZ = 5, Yp = 0.2485 




< 10~^ 


< 10~3 


< 10-3 




7.188 


0.826 


0.355 


AY/AZ = 2, Yp = 0.23 


DL Tau 


15.75 


6.815 


4.181 


UM iau 


19.31 


12.15 


9.790 


AY/AZ = 2, Yp = 0.2485 




14.59 


4.953 


2.782 




16.83 


8.719 


6.555 


AY/AZ = 5, Yp = 0.2485 




9.141 


1.225 


0.515 




6.555 


1.297 


0.675 


AY/AZ = 2, Yp = 0.23 


O Y lau 


5.420 


3.667 


3.256 


BP Tau 


0.180 


0.155 


0.147 


AY/AZ = 2, Yp = 0.2485 




4.675 


3.145 


2.776 




0.138 


0.115 


0.107 


AY/AZ = 5, Yp = 0.2485 




3.058 


2.071 


1.788 




0.054 


0.038 


0.034 


AY/AZ = 2, Yp = 0.23 


GM Aur 


59.83 


9.330 


1.688 


MWC 480 


0.161 


0.162 


0.163 


AY/AZ = 2, Yp = 0.2485 




53.20 


6.241 


0.865 




0.306 


0.309 


0.309 


AY/AZ = 5, Yp = 0.2485 




40.15 


0.817 


0.037 




2.305 


2.325 


2.332 


AY/AZ = 2, Yp = 0.23 


LkCa 15 


3.908 


7.764 


8.856 










AY/AZ = 2, Yp = 0.2485 




5.393 


9.119 


9.496 










AY/AZ = 5, Yp = 0.2485 




10.91 


8.083 


5.268 











M. Gennaro, P. G. Prada Moroni and E. Tognelli 

Table A3. Ba\-('s fac tor \-ahi('s foi' the case of flat prior on the iriass and Gaussiari prior on the niotallicity 



a = 1.20 a = 1.68 a = 1.90 



a = 1.20 a = 1.68 a = 1.90 



AY/AZ -- 
AY/AZ -- 
AY/AZ -- 



: 2, Yp : 
: 2, Yp : 
: 5, Yp : 



0.23 

: 0.2485 
: 0.2485 



RS Cha A 



0.887 
0.752 
0.530 



0.897 
0.752 
0.532 



0.896 
0.764 
0.531 



RS Cha B 



1.058 
0.978 
0.790 



1.089 
1.002 
0.806 



1.114 
1.044 
0.833 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
2, Yp : 
: 5, Yp : 



RXJ A 



0.23 
: 0.2485 
: 0.2485 



0.992 
0.951 
0.873 



1.016 
0.976 
0.905 



1.007 
0.972 
0.904 



RXJ B 



1.353 
1.264 
1.077 



1.071 
0.967 
0.766 



0.975 
0.871 
0.675 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



V1174 Ori A 



V1174 Ori B 



0.23 
: 0.2485 
: 0.2485 



1.103 
1.066 
0.992 



1.034 
0.984 
0.861 



0.991 
0.929 
0.786 



1.313 
1.110 
0.785 



0.999 
0.836 
0.581 



0.924 
0.772 
0.534 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp 
2, Yp 

5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



EK Cep A 



3.545 
2.343 
0.091 



3.568 
2.367 
0.093 



3.570 
2.378 
0.093 



EK Cep B 



2.769 
1.769 
0.783 



0.865 
0.795 
0.699 



0.840 
0.792 
0.710 



AY/AZ : 
AY/AZ : 
AY/AZ 



2, Yp : 

2. Yp 



: 0.23 

0.2185 
0.2f<sr, 



TY CrA A 



TY CrA B 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



ASAS A 



: 0.23 
: 0.2485 
: 0.2485 



0.632 
0.597 
0.544 



0.701 

0.666 
0.610 



0.706 
0.667 
0.604 



ASAS B 



0.840 
0.788 
0.706 



0.680 
0.604 
0.481 



0.581 

0.502 
0.380 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



HD 113449 A 



0.23 
: 0.2485 
: 0.2485 



0.874 
0.986 
1.043 



1.052 
0.948 
0.631 



0.975 
0.827 
0.490 



HD 113449 B 



1.028 
1.325 
1.987 



1.296 
1.636 
2.383 



1.375 
1.731 
2.498 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
: 2, Yp : 
: 5, Yp : 



0.23 
: 0.2485 
: 0.2485 



NTT A 



NTT B 



AY/AZ : 
AY/AZ : 
AY/AZ : 



: 2, Yp : 
2, Yp : 
5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



HD 98800 Ba 



1.499 
1.287 
1.046 



0.676 
0.614 
0.542 



0.560 
0.501 
0.420 



HD 98800 Bb 



1.479 
1.075 
0.843 



0.790 
0.708 
0.595 



0.662 
0.538 
0.425 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



UZ Tau Ea 



UZ Tau Eb 



0.23 
: 0.2485 
: 0.2485 



1.629 
1.446 
1.383 



1.150 
1.006 
0.857 



1.007 
0.845 
0.692 



1.092 
0.687 
0.311 



0.404 
0.252 
0.072 



0.294 
0.161 
0.044 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp 
2, Yp 

5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



DL Tau 



2.058 
1.777 
1.372 



1.052 
0.941 
0.822 



0.892 
0.791 
0.708 



DM Tau 



1.474 
1.248 
0.893 



1.015 
0.853 
0.604 



0.914 
0.767 
0.542 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 
5, Yp : 



: 0.23 

: 0.2485 
: 0.2485 



CY Tau 



1.448 
1.298 
1.131 



1.142 
1.040 
0.829 



1.061 
0.945 
0.731 



BP Tau 



1.105 
1.058 
0.965 



1.078 
1.029 
0.928 



1.068 
1.019 
0.914 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 

2, Yp : 
5, Yp : 



GM Aur 



: 0.23 

: 0.2485 
: 0.2485 



7.278 

5.192 
4.421 



MWC 480 



0.368 

0.218 
0.080 



0.061 
0.027 
0.008 



0.462 

0.313 
0.125 



0.466 

0.316 
0.126 



0.468 

0.317 
0.127 



AY/AZ : 
AY/AZ : 
AY/AZ : 



2, Yp : 
2, Yp : 

5, Yp : 



0.23 
: 0.2485 
: 0.2485 



LkCa 15 



1.243 
1.186 
1.059 



1.048 
0.965 
0.791 



0.959 
0.869 
0.690 



Testing pre-main sequence models: the power of a Bayesian approach. 

Table A4. BaycK fac tor \-ahi('s foi' the case with Gaiissian prior on liotli the mass and the inotallicity 







a = 1.20 


a = 1.68 


a = 1.90 




a = 1.20 


a = 1.68 


a = 1.90 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


Ho Una A 


17.40 
10.23 
0.003 


17.78 
10.24 
0.003 


17.51 
10.23 
0.003 


tib Ona D 


26.67 
10.36 
0.003 


27.27 
10.60 
0.003 


27.57 
10.85 
0.003 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


RXJ A 


10.87 
7.781 
3.366 


18.35 
18.34 
14.13 


17.02 
18.31 
15.14 


RXJ B 


15.62 
15.98 
12.00 


14.00 
10.44 
3.298 


11.73 
7.579 
1.666 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


V1174 Un A 


2.484 
6.836 
22.02 


12.76 
18.37 
13.08 


16.35 
18.09 
7.076 


V1174 (Jri B 


0.794 
0.217 
0.007 


0.091 
0.016 
< 10-3 


0.043 
0.006 
< 10-3 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


EK Cep A 


37.95 
8.021 
0.001 


38.17 
8.004 
0.001 


38.06 
8.010 
0.001 


EK Cep B 


24.89 
16.38 
4.962 


21.46 
28.07 
16.00 


9.717 
17.43 
22.97 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2. Yp = 0.2185 
AY/AZ = 5p = 0.21,sr, 


1 Y Or A A 


- 


- 


- 


1 Y UrA D 


- 


- 


- 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


A O A O A 

AbAb A 


3.527 
10.47 
26.83 


1.714 

4.188 
10.93 


3.897 
6.801 
11.13 


AoAo B 


13.81 
13.61 
7.133 


4.101 
1.690 
0.056 


1.250 
0.340 
0.003 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


tiu 113449 A 


8.426 
7.398 
3.911 


9.087 
6.058 
1.671 


7.972 
4.900 
1.155 


T TTA ~i ~i ^ A A r\ T~> 

HD 113449 B 


4.084 
4.029 
3.304 


5.100 
4.877 
3.765 


5.393 
5.124 
3.882 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


JN i i A 


- 


- 


- 


■\Trprp T~) 
JN 1 i B 


- 


- 


- 


AY/AZ = 2, Yp = 0.23 
AY/AZ = 2, Yp = 0.2485 
AY/AZ = 5, Yp = 0.2485 


HD 98800 Ba 


4.598 
5.117 
5.514 


4.732 
3.838 
2.523 


3.154 
2.254 
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Figure Bl. RS Cha components mass and age distribution functions as obtained from the logg — logTcfi diagram using the standard 
set of models. Upper panels: Marginalized mass distributions. Lower panels: Marginalized age distributions. Left panels: Marginalization 
using a flat mass prior. Right panels: Marginalization using a Gaussian mass prior. All panels: In dot-dashed-blue the primary component, 
in dashed-red the secondary. Full symbols indicate the mode of the distributions, the bars mark the 68% confidence interval. The quoted 
numbers represent the relative precision of the mass or age estimates. In the upper panels the empty symbols and related error bars 
indicate the dynamical masses and their measurement errors. In the lower panels the solid-black line represents the system age distribution. 
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Figure B2. RXJ 0529.4+0041 components mass and age distribution functions from comparison with the standard set of models (see 
Fig. lBll for a description). 
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Figure B3. V 1174 Ori components mass and age distribution functions from comparison with the standard set of models (see Fig. IBll 
for a description). 
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Figure B4. V 1174 Ori components mass and age distribution functions from comparison with the standard set of models (see Fig. IBll 
for a description). In this case the temperature of the primary has been artificially raised by 100 K and the temperature of the secondary 
has been raised accordingly in order to keep the temperature ratio constant. 
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Figure B5. EK Cep components mass and age distribution functions from comparison with the standard set of models (see Fig. lBll for 
a description). 
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Figure B6. EK Cep components mass and age distribution functions from comparison with the set of models with a = 1.20, Yp = 0.23 
and AY/ AZ = 2 (see F ie:.IBl|for a description) . The stellar Z values used for this comparison have been calculated using the observed 
[Fc/H] and {Z/X)q bv lCrevesse fc Sauvallfiggsl . 
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Figure B7. TYCrA components mass and age distribution functions from comparison with the standard set of models (see Fie. lBll for 
a description). 
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Figure B8. TY CrA components mass and age distribution functions from comparison with the set of models with a = 1.20, Yp = 0.23 
and AY/AZ = 2 (see Fig. IB II for a description). 




Figure B9. ASAS J052821+0338.5 components mass and age distribution functions from comparison with the standard set of models 
(see Fig. IBll for a description). 
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Figure BIO. ASAS J052821+0338.5 components mass and age distribution functions from comparison with the standard set of models. 
In this case we used the data from Table 1 of iStempels et al..i200S when no spots are included in the light-curve modeling (see Fig. IBll 
for a description). 




Figure Bll. HD 113449 components mass and age distribution functions from comparison with the standard set of models (see Fig. IBll 
for a description). 
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Figure B12. HD 113449 components mass and age distribution functions from comparison with the set of models with a = 1.20, Yp 
0.23 and AY/AZ = 2 (see Fig. iBllfor a description). T l 
observed [Fe/H] and {Z/X)q bv lcrevesse fc Sauvallfiggsl . 
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Figure B13. NTT 045251+3016 components mass and age distribution functions from comparison with the standard set of models (see 
Fig. IB II for a description). 



38 M. Gennaro, P. G. Prada Moroni and E. Tognelli 




0.2 0.4 



0.6 0.8 1.0 



1.2 1.4 1.6 






0.222 - 
0.250 - 
0.642 i 













0.199 - 






0.217 - 






0.506 ; 






















' ^ \ 




: I 


1 ^ \ 




■ ' 







10 15 20 

r [Myr] 



10 15 20 

r [Myr] 



Figure B14. NTT 045251+3016 components mass and age distribution functions from comparison with the set of models with 
1.20, yp = 0.23 and AY/AZ = 2 (see Fig.[BT]for a description). 
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Figure B15. HD 98800 B components mass and age distribution functions from comparison with the standard set of models (see Fig 
for a description). 



